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This paper considers a GI/GI/1 processor sharing queue in which 
jobs have soft deadlines. At each point in time, the collection of resid- 
ual service times and deadlines is modeled using a random counting 
measure on the right half-plane. The limit of this measure valued pro- 
cess is obtained under diffusion scaling and heavy traffic conditions 
and is characterized as a deterministic function of the limiting queue 
length process. As special cases, one obtains diffusion approximations 
for the lead time profile and the profile of times in queue. One also 
obtains a snapshot principle for sojourn times. 

1. Introduction. Congestion and delay is typical in heavily loaded queue- 
ing systems and may vary significantly according to the service discipline 
being used. There is much interest in identifying service disciplines that 
minimize delay and this has led to a rich literature on performance analy- 
sis. Classical measures of delay include workload, queue length and sojourn 
time. In many circumstances, these performance measures help assess a ser- 
vice discipline's effectiveness at minimizing the time jobs spend in the sys- 
tem (their sojourn times). However, these measures may be inadequate for 
systems with more specific timing requirements. For example, in systems 
which have heterogeneous job deadlines, minimizing sojourn times may not 
be the best strategy for minimizing missed deadlines. Recent work has in- 
vestigated various queueing models that allow for a more general notion of 
delay [4, 6, 13, 14, 15, 16, 18, 26]. In these models, jobs have individual 
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random deadlines and a service discipline's performance depends on how 
effectively it meets them. 

Such real-time queueing models arise naturally in applications such as 
manufacturing systems, voice and video communication systems and avion- 
ics and automotive control systems. The associated performance measures 
are high-dimensional because individual timing information must be tracked 
to determine if deadlines are being met. Exact analysis is usually intractable, 
but heavy traffic approximations often exist. This was verified for the EDF 
(earliest deadline first) service discipline in [6] for the GI/GI/1 queue, in 
[15] for the multiclass queue with HLPS (head-of-the-line processor sharing) 
across classes, in [26] for feed-forward networks and in [16] for multiclass 
acyclic networks. Corresponding results for the FIFO (first in first out) dis- 
cipline were also given in [15, 26]. The accuracy of these approximations in 
the single server case was investigated in [13, 14]. 

Processor sharing is a widely used idealization of the round-robin and 
time-sharing protocols used in computer and communication systems. Due 
to its practical importance, it is natural to investigate analogues of the 
results in [6, 15, 16, 26] for this service discipline. This is the subject of the 
present paper, which provides a framework for evaluating the heavy traffic 
performance of the GI/GI/1 processor sharing queue with respect to a fairly 
general structure of deadlines. 

A new feature of this framework is that job deadlines are allowed to be 
correlated with service times; previous work has assumed independence of 
the two. This makes it possible, for example, to model scenarios in which 
large jobs have longer deadlines than small jobs. Although this adds a level 
of realism, it also requires some additional technical machinery. The state of 
the processor sharing queue with deadlines will be tracked using a measure 
valued process in the right half-plane. This idea builds on previous work 
on the GI/GI/1 processor sharing queue [8, 9, 10, 20, 21]. The setup re- 
quires detailed information about how arriving jobs affect the system state. 
Consequently, the primitive processes must be considered jointly as a mea- 
sure valued arrival process. This creates some difficulty. The usual estimates 
needed to control scaling limits of an arrival process derive from functional 
weak laws of large numbers and functional central limit theorems. Analo- 
gous estimates for the arrival process used here require stronger results: the 
paper will draw on the theory of empirical processes to obtain functional 
Glivenko-Cantelli-type estimates for the measure valued arrival process (see 
Section 4.2). 

Although dealing with measure valued processes in the half-plane requires 
some theoretical overhead, the additional work pays off. Section 3 provides 
some examples of the types of computation the more general setup allows. 
In particular, Section 3.2 uses it to compute a processor sharing analogue of 
Reiman's snapshot principle [18, 22] and Section 3.3 gives an example of how 
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one might set realistic deadlines under the assumption of linear dependence 
between service times and deadlines. Furthermore, the present analysis for 
processor sharing is just an example of how empirical process theory may 
be applied to measure valued state descriptors. It is likely that the approach 
can be adapted to the study of other models. 

The model considered here consists of an infinite capacity buffer, to which 
jobs arrive according to a delayed renewal process. Jobs arrive with service 
time requirements determined by a sequence of independent, identically dis- 
tributed positive random variables. A single server processes buffered jobs 
at unit rate according to the processor sharing discipline: it works simulta- 
neously on all jobs in the buffer, providing an equal fraction of its capacity 
to each. That is, if there are Z(t) > 1 jobs in the buffer at time t, then each 
job is receiving service at the instantaneous rate 1/Z(t). When the server 
has fulfilled a given job's service time requirement, the job exits the system. 
In addition to its service time, each job arrives at the buffer with a deadline 
given by the job's arrival time plus a real valued random variable called the 
initial lead time. The sequence of initial lead times is independent and iden- 
tically distributed, but a job's initial lead time may be correlated with its 
service time. Job deadlines are soft, meaning that jobs remain in the system 
until served to completion; a job still in the system when its deadline expires 
is called late. 

Since deadlines are soft and the server does not take them into account, 
the processor sharing discipline is unaltered by the inclusion of timing infor- 
mation in the model. In particular, the workload and queue length processes 
are identical to those of the classical processor sharing model. Thus, the re- 
sult described below should be interpreted as a performance analysis of clas- 
sical processor sharing, where the performance is measured with respect to 
the aforementioned deadline structure. In particular, the result answers the 
following question. Given a certain structure of deadlines, how well does the 
processor sharing discipline perform in meeting them? The setup used al- 
lows a variety of functionals to be considered in answering this question. For 
example, the lead time profile, time-in-queue profile and empirical sojourn 
time distributions can all be described using the main result (see Section 3). 
In this regard, the present analysis may be used to compare the effectiveness 
of processor sharing with that of other disciplines under the same deadline 
structure. As mentioned above, performance analyses of this type exist for 
certain FIFO models (another discipline that ignores the deadlines) as well 
as certain EDF models (a discipline which takes them into account). 

Note that the system studied here is different from a GI/GI/1 processor 
sharing queue with firm deadlines, in which jobs exit early if their deadline 
expires. The latter system exhibits very different behavior and is the subject 
of forthcoming work. For more on processor sharing queues, see [25] for a 
survey up to 1987 and [10] for a discussion of more recent work. 
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The main result of this paper concerns a measure valued process that 
tracks the state of the processor sharing queue with soft deadlines. At time 
t > 0, this state includes both the residual service time and the lead time of 
each buffered job: if I(i) denotes the set of jobs in the buffer at time t, then 
the residual service time Vi{t) is the remaining service requirement of job 
i E I(t) and the lead time li(t) is the deadline of job i minus the current time 
t. If li{t) is non-negative, it represents the remaining time until the deadline 
of job i expires; if it is negative, its absolute value is the time overdue. Let 
M be the set of finite, non-negative Borel measures on the right half-plane 
H + = [0,oo) x (—00,00) and let Sr x y -\ denote the Dirac point measure at 
(x,y) E H+. For each t > 0, the state descriptor Z{t) is the random element 
of M given by 

iei(t) 

Note that projection of Z(-) onto the first coordinate yields the state de- 
scriptor used in [8, 9, 10, 20, 21]. For each r in a sequence 1Z of positive 
real numbers tending to infinity, define a diffusion scaled version of the state 
descriptor by 

Z r {t)(B x C) = -Z(r 2 t)(B x rC), Bc[0,oo), Cc (-00,00). 
r 

Let a be the limiting arrival rate (as r — > 00) for jobs entering the system 
and let $ E M be the limiting joint distribution (as r — ► 00) of service times 
and initial lead times (the precise form of these assumptions is given in 
Section 2.3 below). For z > 0, let i?g E M be the measure defined by 

$l([x, 00) x [y, 00)) = a "d([x + uz~ l , 00) x [y + u, 00)) du 

Jo 

for all (x,y) E H + and let $® = 0, where denotes the zero measure in M. 
The main result of this paper states that under mild conditions, including 
standard heavy traffic assumptions, the sequence of diffusion scaled state 
descriptors {Z r (-)} converges in distribution as r — ► 00 to the measure valued 

process C w , where Z*(-) is a reflected Brownian motion on [0, 00). The 
process Z*(-) is the weak limit of the diffusion scaled queue length process 
obtained in [9]. Note that for each z > 0, the measure can be computed 
in a straightforward way, either numerically, via simulations, or, in some 

cases, explicitly. Thus, the process v e is a tractable approximation to 
the dynamics of the processor sharing queue with soft deadlines. 

The proof of this result proceeds in several stages and is summarized 
informally as follows. For each t > 0, define a fluid scaled version of the 
state descriptor by 

Z r (t)(B xC) = -Z(rt){B x rC), B C [0, 00), C C (-00, 00), 
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and for m = 0, 1, 2, . . . , define shifted processes Z r ' m (t) = Z r (m + t). Fix 
times T, L > 1 . It is first proved that 7£-indexed sequences of sample paths 
of the shifted processes {Z r,m {t) :r G TZ,m < [rT\,t G [0,L)} are precom- 
pact with high probability. Next, limit points of these sequences are charac- 
terized as local fluid limits, which are continuous measure valued functions 
£ : [0, L) — > M. In particular, it is shown that these local fluid limits uniformly 
approximate the sample paths of Z r ' m (-) on [0,L) with asymptotically high 
probability. Moreover, for sufficiently large L and t G [L — 1,L), every local 
fluid limit £(•) is approximately in steady state, which is given by 

(1.1) C(*)«^ (t) , 

where z(t) = ((t)(M + ) is the total mass of £(t). 

For each r G 1Z, the interval [0, rT] is covered by the overlapping time in- 
tervals [to, 771+ L), where m = 0, . . . , [rT] . Thus, for any t G [0, rT] , there ex- 
ist to < [rT\ and s G [0,L) such that Z r (t) = iT> m (s). Since Z r (t) = Z r (rt), 
the sample paths of {Z r ' m (-) :m < L r ^J} ° n the time interval [0, L) deter- 
mine the sample paths of Z r (-) on [0,T]. Consequently, (1.1) translates to 
the sample paths of Z r (-) on [(L — l)r _1 ,T]. That is, as r — > oo, 

(1.2) Zr{t)^^\ te[(L — l)r _1 ,T], 

where -Z r (-) = ^ r (-)(IH + ) is the diffusion scaled queue length process. With 
some extra work, one obtains (1.2) for all t G [0, T]. 

In particular, the diffusion scaled state descriptor Z r (-) can be asymp- 
totically recovered from the one-dimensional diffusion scaled queue length 
process Z r (-) by the lifting map A^:zh>^. This phenomenon is known 
as state space collapse. The process Z r {-) converges in distribution to a 
reflected Brownian motion Z*(-) by Corollary 2.4 of [9]. Applying the con- 
tinuous mapping theorem to (1.2) and A# completes the proof. 

The program outlined above for proving state space collapse of a dif- 
fusion scaled process using overlapping sections of a fluid scaled process 
is motivated by methodology developed in [3, 24] for open multiclass net- 
works with HL (head-of-the-line) service disciplines. This methodology was 
adapted to the measure valued state descriptor of a processor sharing queue 
in [9] and has been further adapted here. The limiting diffusion process 

A$Z*[-) = i? e takes values in a submanifold {t?e : z > 0} of M, known as 
the invariant manifold. Note that this manifold is not a linear space. This 
is an interesting contrast to other results exhibiting state space collapse, 
including [3, 5, 9, 24], in which the invariant manifold is frequently a linear 
space. 

The paper is organized as follows. Section 2 gives a precise description of 
the model, different forms of scaling and various asymptotic assumptions; 
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it also contains the statement of the main theorem. Section 3 illustrates 
several applications of the theorem. Section 4 contains the proof of tightness 
of the shifted fluid scaled state descriptors. Section 5 investigates properties 
of local fluid limits and Section 6 contains the proof of state space collapse, 
leading quickly to the proof of the main theorem. 

1.1. Notation. The following notation will be used throughout. Let N = 
{1, 2, . . .} and let Z and R denote the integers and real numbers, respectively. 
Let R + = [0, oo) and denote the right half-plane R + x R by H+. For a,b G 
R, write a V b for the maximum, a A b for the minimum, a + and a - for 
the positive and negative parts, [a\ f° r the integer part and [a] for the 
smallest integer n> a. Denote the indicator of a set B C H+ by 1b and let 
B e = {w £ H + : mi z( zB \\w — z\\ < e}. For w G H + , let B + w = {z + w: z£ B}. 
Denote by C&(<S) the space of continuous bounded real valued functions on 
a topological space S. 

Recall that M is the set of nonnegative finite Borel measures on W+, 
with zero measure denoted by 0. For £ G M and a £-integrable function 
g : H_|_ — > R, define (g, = J H+ g dQ. The space M is endowed with the weak 

topology: Cn C in M if and only if {g,(n) -> (g,Q for all g G C b (H + ). 
With this topology, M is a Polish space [19]. It will be convenient to use the 
following metric: for £, £ G M, let 

d[£, C] = inf{e > : £(B) < ((B e ) + e and 

(1.3) 

C(B) < i{B £ ) + e, for all closed B C H+}. 

It is straightforward to verify that d[-, •] is a complete metric on M inducing 
the weak topology. 

All stochastic processes are assumed to be right continuous with finite 
left limits. Let D([0, oo), M) denote the space of right continuous, left lim- 
ited paths C : [0, oo) — ► M, endowed with the Skorohod J\ -topology. For two 
random objects X, Y with the same distribution fj,, write X ~ Y ~ fi. Write 
X n Y when X n converges in distribution to Y as n — > oo. 

2. A processor sharing queue with soft deadlines. This section gives a 
precise description of the model under consideration, specifies assumptions 
and states the main result. A formal definition of the processor sharing 
queue was given in [10] and the reader is referred there for a more detailed 
description. The present model is a generalization of the one considered 
in [10], so parts of the definition are restated here. This is the subject of 
Section 2.1 below, which introduces a sequence of processor sharing models 
incorporating soft deadlines, along with the associated notation. Section 2.2 
describes the scaling and time shifts to be applied to this sequence and 
Section 2.3 specifies asymptotic assumptions. The main result appears in 
Section 2.4. 
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2.1. Sequence of models. Let 1Z C (0, oo) be a sequence that increases to 
infinity. Suppose that, for each r G 1Z, there is a stochastic model consist- 
ing of the following: a processor sharing server that processes jobs at unit 
rate from an infinite capacity buffer, a collection of stochastic primitives 
(E r (-), {v\, describing the arrival times, service times and initial lead 

times of jobs arriving to the buffer and a random initial condition speci- 
fying the state of the system at time 0. The time evolution of the system 
state is described by a collection of performance processes, defined in terms 
of the primitives and initial condition through a set of descriptive equa- 
tions. The random objects in each model are defined on a probability space 
(f2 r , J^ r , P r ), with expectation on this space denoted E r . 

The stochastic primitives consist of an exogenous arrival process E r {-) 
and a sequence of service times and initial lead times {v\, V i } c ?l l . The arrival 
process E r (-) is a rate a r delayed renewal process associated with a sequence 
{u^^i of finite nonnegative interarrival times. For t > 0, E r {t) represents 
the number of jobs that have arrived at the buffer during the time interval 
(0,t]. The ith job to arrive after time is called job i; jobs already in the 
buffer at time will be called initial jobs. The quantity u\ is the arrival 
time of the first job and for i>2, u\ is the elapsed time between the arrival 
of job i — 1 and job i. Thus, job i arrives at time U[ = Y?j=i uT j f° r £ > 1- 
Define E/jJ = 0. Then for t > 0, 

E r (t)=sup{i>0:Ul' <t}. 

Assume that {u r i }'^ l is a sequence of independent random variables and that 
{u[}^ 2 are independent and identically distributed with mean l/a r G (0, oo) 
and standard deviation a r < oo. The first element u[ of the sequence is 
assumed to be strictly positive with finite mean. 

For each i > 1, the service time v\ represents the amount of processing 
time that job i requires from the server; the initial lead time represents 
the maximum amount of time that job i can be in the buffer without being 
late. Since this job arrives at time UJ, it will be late after time U[ + 1\. Late 
jobs remain in the buffer until completing service. Consequently U[ + l\ 
is a soft deadline for job i. Assume that {v\} are strictly positive random 
variables and that is a sequence of independent and identically 

distributed random vectors with common joint distribution given by a Borel 
probability measure # r on H+. Assume that v\ has mean 1/ (3 r G (0, oo) and 
standard deviation b r < oo and that {v^,l^}°2 =1 is independent of the arrival 
process E r {-). 

It will be convenient to express the primitives as a measure valued arrival 
process. For each t > 0, let £ r (t) G M be the random measure 

ET(t) 
i=l 
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where S( x ,y) denotes the Dirac measure at (x,y) € H + . 

The initial condition specifies Z r (0), the number of initial jobs present in 
the buffer at time zero, as well as the service time requirements and initial 
lead times of these initial jobs. Assume that Z r (0) is a nonnegative, integer 
valued random variable. The service times and initial lead times for initial 
jobs are the first Z r (0) elements of a sequence {vj , lj}'^ =1 of random vec- 
tors, where {vj} are strictly positive. Assume that Z r (0) and {uJ,/J}^ =1 are 
independent of {m[}^ 2 an d { v l Ji}i^=i\ they do not need to be independent 
of each other. A convenient way to express the initial condition is to define 
an initial random measure Z r (0) € M by 

z r (o) 

Z r (0)= J2 S^jry 
3=1 

Henceforth, Z r (0) will be used as the initial condition. Let x ; H+ ~ > ^+ 
denote the projection (x,y) t— > x. Assume that Z r (Q) satisfies 

(2.1) W[(l,Z r (0))]<oo, 

(2.2) E r [{ X ,Z r (0))]<oo. 

Note that since (l,Z r (0)) = Z r (0) and ( x , Z r {0)} = Y^fL? 6J, assumptions 
(2.1) and (2.2) mean that the expected initial queue length and expected 
initial workload are finite. 

The performance processes describe the system's behavior; they are de- 
fined in terms of the primitives through a set of equations that embody the 
processor sharing discipline and the expiration of deadlines. For each t > 0, 
let S r (t) denote the cumulative service per job provided by the server up to 
time t. Thus, if job i < E r {t) arrives at time U[ < t, the cumulative amount 
of processing time job i receives by time t is equal to v\ A (S r (t) — S r (U[)). 
Similarly, the cumulative amount of processing time an initial job j < Z r (0) 
receives by time t is equal to Vj A S r (t). Define the residual service time at 
time t of job i (and initial job j) by 

(2.3) <(i) = « - S r (t) + S r (Un) + , i = l,...,E r (t), 

(2.4) v r j (t) = (v r j -S r (t)) + , j = l,...,Z r (0). 

A job's residual service time represents the current remaining amount of 
processing time required to fulfill its service time requirement. If a job's 
residual service time is zero, it has completed service and departed the buffer. 
Define the lead time at time t of job i (and of initial job j) by 

(2.5) m = U r + J r_ t ^ i = li ... iE r(t), 

(2.6) l]{t) = V 3 -t, J 1 Z'iO). 
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The lead time of a job currently in the buffer is the time remaining until its 
deadline passes; jobs with negative lead times are currently late. 

The state descriptor Z r (-) : [0, oo) — > M is a measure valued performance 
process that describes the time evolution of the system state: at each time t > 
0, the random measure Z r (t) has one unit of mass located at (vj(t),lj(t)) 6 
H + for each initial job j < Z r (0) still in the buffer at time t and one unit of 
mass at /£(£)) £ H+ for each job i < E r (t) still in the buffer at time 

t. For (x,y) G H+, let be the Dirac measure at (x,y) if x > 0, with 

^(o y) = ^' Since a job is still in the buffer if and only if its current residual 
service time is positive, we have 

Z r (0) E r (t) 
j=l i=l 

Let ^ r (t) denote the number of jobs in the buffer, or queue length, at time 
t > 0. Then Z r (t) = (l,i? r (i)). Thus, under the processor sharing discipline, 
any job present in the buffer at time t receives service at the instantaneous 
rate (l,Z r (t))~ 1 . Note that if a job is present in the buffer at time t, then 
(l,Z r (t)) / 0. Let <p(x) = l/x for x E (0,oo), with ip(0) = 0. Then the cu- 
mulative service per job up to time t can be written 

(2.8) S r (t)= f\{(l,Z r {s)))ds. 

Jo 

It will be convenient to let t = S r (t) — S r (s) for t > s > 0. 

Given the primitives (E r (-), {v\, l\} £2^) and the initial condition Z r (0), 
the equations (2.3)-(2.8) determine the residual service times, lead times, 
the cumulative service per job process S r (-) and the state descriptor Z r (-). 
This fact is not difficult, although somewhat tedious, to show. 

Let W r {t) denote the workload in the buffer at time t > 0. This is the 
amount of time the server would have to work to complete the remaining 
service time requirements of all jobs in the buffer at time t, assuming no new 
arrivals take place. Since this equals the sum of the residual service times of 
all jobs present in the buffer at time t, we have 

W r (t) = ( X ,Z r (t)). 

Note that the sequence of processor sharing models introduced here is a 
generalization of the sequence of models studied in [9, 10]. In particular, the 
measure valued process //(•) = Z r {-) o \~ l is determined by the primitives 
(£ r (•),{<},=!) and the initial condition (Z r {ti),{v^}f =1 ) and is the state 
descriptor that was studied in [9, 10]. 
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2.2. Scaling. This paper concerns the asymptotic behavior, under heavy 
traffic conditions, of the 7?.-indexed sequence of models introduced in Sec- 
tion 2.1. To obtain useful limits as r — > oo, various objects comprising the 
rth model must be appropriately scaled. 

In Section 2.3 below, it is assumed that, as r — > oo, the processor sharing 
models become heavily loaded at a rate governed by r. This implies that a 
job in the rth model remains in the buffer for a time that is of the order r 
multiplied by the job's service time requirement. Thus, if initial lead times 
were to remain of the same order as service times, lead times l\{t), for large 
t (of order r 2 ), would tend to — oo as r — > oo. To obtain nontrivial scaling 
limits for the lead times, initial lead times in the rth model are assumed 
to be of order r and will be scaled, along with lead times by r . For 
each r G 1Z, let t? r G M be the probability measure satisfying 

(2.9) ti r {B xC)=tf r {B xrC) 

for all Borel sets B C R+ and C C R. 

The asymptotic heavy traffic behavior of the state descriptor will be exam- 
ined on diffusion scale. For each r G 1Z, the diffusion scaled state descriptor 
is defined, for t > 0, as the random measure Z r (t) G M satisfying 

Z r (t)(B xC) = -Z r (r 2 t)(B x rC) 
r 

for all Borel sets B C M+ and CcK. Note that this definition also scales 
lead times by r . 

Diffusion-scaled versions of the workload and queue length processes are 
also needed. For t > 0, define W T {t) = r^W^rH) and Z r (t) = r- x Z r (rH); 
note that W r (t) = (x,2 r (t)) and Z r {t) = (l,Z r (t)). Since the workload 
process of a single server queue is the same for all work conserving ser- 
vice disciplines, including processor sharing, the heavy traffic behavior of 
W r (-) is described by the well-known result for FIFO queues [11]. Since 
Z r (-) = (l,Z r (-)) = (l,Z r (-) o x _1 ), the queue length process of the present 
model is the same as for the processor sharing queue without deadlines; a 
heavy traffic analysis of Z r {-) appears in [9]. The existing results for W r (-) 
and Z r (-) will be used in the analysis of Z r {-). 

Results for the diffusion scaled state descriptor Z r (-) will be derived from 
results for the fluid scaled state descriptor, defined for t > as the random 
measure Z r {t) G M satisfying 

(2.10) Z r (t)(B xC) = -Z r (rt)(B xrC) 

r 

for all Borel sets BcK + and CcR. 

The relationship Z r (t) = Z r {rt) will be essential to bootstrapping results 
from fluid scale up to diffusion scale. The process Z r {-) is considered over a 
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fixed finite time interval [0, T] for T > 1. This corresponds to looking at the 
fluid scaled process Z r (-) over [0,rT]. As r — > oo, the asymptotic behavior 
of Z r {-) over finite time intervals is analyzed using modified versions of 
the techniques developed in [9, 10]. Since [0,rT] grows without bound as 
r — > oo, it is necessary to piece together many (order r) overlapping sections 
of Z r (-), each defined on a finite time interval of fixed length L > 1. This 
strategy is analogous to the one used in [3]. These overlapping sections of 
Z r (-) are called the shifted fluid scaled state descriptors. For each t > and 
m G {0, 1, . . .}, define 

(2.11) Z r ' m {t) = Z r {m + t). 

Then for each r G 1Z, the time interval [0,rT] is covered by the \rT\ + 1 
overlapping time intervals [m, m + L) for m = 0, . . . , \rT\ . Observe that for 
each t G [0,rT], there exists (at least one) m < [rT\ and s £ [0,L] such that 

Z r {t) = Z r ' m {s). 

Analysis of the processes Z r,m {-) will involve fluid scaled and shifted fluid 
scaled versions of many of the processes introduced so far. For all r € 1Z, 
me {0,1,...}, t>0 and i = l,...,E r (rt), define 

(2.12) E r (t) 

(2.13) S r (t) 

z r (t) 

(2.14) vl{t) 

(2.15) l\{t) 

For r G TZ and t > 0, define the fluid scaled measure valued arrival process 
by 

(2-16) C r {t) = - £ V.irr-1)- 

i=i 

The fluid scaled processes S r (-) and >C r (-) will play particularly important 
roles; it will be convenient to have notation for their increments. Define the 
following fluid scaled and shifted fluid scaled increments: for all r G TZ, all 
t > s > and all m G {0, 1, . . .}, let 

(2.17) Sl t = S r (t)-S r (s), S r J 

(2.18) Cl t = C r (t)-C r (s), 



= -E r (rt), E r ' m (t) = E r (m + t), 
r 

= S r {rt), S r ' m (t) = S r (m + t), 

= -Z r {rt), Z r ' m {t) = Z r {m + t), 
r 

= vl{rt), ^ m {t)=vl(t + m), 

= -ll{rt), f> m {t)=ll(t + m). 



S 
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Note that by (2.17) and (2.13), 

rr{m+t) rt 

S r s f= / V ((l,Z r (u)))du= / V ((l,Z r > m (u)))du. 

Jr{m+s) Js 

2.3. Asymptotic assumptions. This section imposes asymptotic assump- 
tions on the sequence of models introduced in Section 2.1. This is the setting 
in which the main theorem is proved. 

Let a, a and p be fixed positive constants and let 7 € K. Let 1? be a 
probability measure on H+ satisfying 

(2.19) i?({0}xl)=0, 

(2.20) <x 4+ V)<oo, 

(2.21) ( X ,#)=a~ 1 . 

Then 6= ((x 2 ,i9) - (x,$) 2 ) 1/2 is finite. For the sequence of arrival processes, 
assume that as r — ► 00, 

(2.22) (a r ,a r ) ^(a,a), 

(2.23) E r [<]/r^0, 

(2.24) limsu P E r [(^) 2+p ] < 00. 

r— >oo 

For the sequence of service times and initial lead times, assume that as 
r — > 00, 

(2.25) tf r -^tf, 

(2.26) (P r ,b r ) - (a, 6), 

(2.27) limsup(x 4+ V r ) < 00. 

r— +00 

Define the traffic intensity parameter for the rth system by p r = a r / (3 r . 
Assumptions (2.22) and (2.26) imply that p r — > 1, that is, the sequence of 
models approaches heavy traffic. Assume, further, that as r — > 00, 



r n 



(2.28) r(l - p 

Assumption (2.28) represents the usual requirement that the sequence of 
systems approaches heavy traffic at a particular rate; the constant 7 ap- 
pears as the drift coefficient of the limiting workload process. Assumption 
(2.23) implies that the initial residual interarrival time vanishes on diffusion 
scale. Assumption (2.25) specifies a weak limit for the joint distributions of 
service times and initial lead times when scaled according to (2.9). Assump- 
tions (2.24) and (2.27) imply Lindeberg-type conditions which, along with 
(2.22), (2.23), (2.26) and (2.28), guarantee functional central limit theorems 
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for the triangular arrays {u\\ z = 1,2,.. .} r g7?. and i = 1, 2, . . .} r( =7j and a 
functional Donsker theorem for the measure valued arrival processes C r {-). 
The central limit theorems ultimately imply that the diffusion scaled work- 
load and queue length processes converge in distribution to reflected Brown- 
ian motions (see Propositions 4.1 and 4.2). The functional Donsker theorem 
implies a Glivenko-Cantelli-type estimate for C r (-). Assumption (2.27) is 
two moments stronger than is normally necessary for a functional central 
limit theorem. The additional restriction is required by Corollary 2.4 in [9], 
which will be applied in the present paper; see [9] for further discussion. 

It remains to impose asymptotic assumptions on the diffusion scaled initial 
condition Z r (0). They are stated in terms of Z r (0) = Z r (0) since they will 
be used in that form. The following definition is needed. 

Definition 2.1 (Invariant manifold). For each z > 0, let i^eMbe the 
unique measure satisfying 

POO 

(2.29) (l[a;,oo)x[i/,oo)i#e) = a \ (^[x+uz- 1 ,oo)x [y+u,oo) > du 

J 

for all x £ and j/ER. Let ft® = and define 

M* = {^eM:z>0}. 

By (2.21), (2.29) and a change of variables, (l,#e) = z f° r an z >0- For 
z > 0, the probability measure z~ lr &l can be thought of as the excess life- 
time distribution of $ in direction (z -1 , 1), a generalization of the notion of 
excess lifetime distribution to the half-plane H+. Following usage in [3], we 
call the one-parameter family of measures C M the invariant manifold 
associated with Note that the invariant manifold is not a linear space. 
Let 6 £ M be a random measure such that 

(2.30) G G M,? a.s., 

(2.31) E[(l,9)]<oo. 

For the sequence of scaled initial measures Z r (0) =Z r (0), assume that as 

r — > oo, 

(2.32) (^(o),(x,^(o)),(x 1+p ,^(o)))^(e,(x,e),( x 1+p ,e)). 

The second component of (2.32) implies that the scaled initial workload 
converges in distribution, that is, W r (0) =>• (x, ©")■ 

The assumptions of this section are now summarized to simplify the state- 
ment of results that follow. 
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There is a sequence of processor sharing models, as 
defined in Section 2.1; there exist positive constants 
(A) a, a, p, a real number 7 £ M, a probability measure 

on H + and a random measure G M such that 
(2.19)-(2.32) hold. 

2.4. Limit theorem. Assume (A) and let Z*(-) be a reflected Brownian 
motion on R + with drift — 27a ( 1 + a 2 b 2 ) ~ 1 , variance 4a 3 (a 2 + b 2 ) ( 1 + a 2 b 2 ) ~ 2 
and initial value Z*(0) that is equal in distribution to (1,0). Define the 
measure valued process 



Theorem 2.2. As r — > 00, i/ie sequence of diffusion scaled state descrip- 
tors {Z r (-)} converges in distribution to the measure valued process Z*{-). 

Sections 4-6 are devoted to proving Theorem 2.2. Section 4 establishes 
a tightness property for the families of shifted fluid scaled state descriptors 
{Z r ' m { ) :m < [rT\}. This yields the existence of limit points, called local 
fluid limits, which are characterized in Section 5 as solutions of a certain 
integral equation. In Section 6, the steady state behavior of these local fluid 
limits, combined with the overlapping nature of the shifted fluid scaled state 
descriptors {Z r ' m (-) :m < [rT\}, reveals that {Z r (-)} is asymptotically close 
to the invariant manifold M$. Combined with an existing limit theorem for 
the diffusion scaled queue length processes {Z r (-)}, this establishes Theo- 
rem 2.2. The next section illustrates the applicability of Theorem 2.2 by 
discussing several useful computations. 

3. Special cases and applications. A central motivation of this paper is 
to study the heavy traffic performance of the processor sharing discipline 
in relation to timing requirements. A useful performance measure in this 
context is called the lead time profile (see [6, 13, 14, 15, 16, 18, 26]). If 
7r:H + — ► M is the projection (x,y) 1— ► y, then the diffusion scaled lead time 
profile is 

Z r {-)o7T-\ 

It is a random finite Borel measure on R that describes how well the service 
discipline is meeting the soft deadlines of jobs. For example, the (scaled) 
number of jobs that are currently late at time t is (lf OOj0 i,Z r (t) o ir~ 1 ). 

This section discusses a few special cases for which the lead time profile is 
given by a convenient formula. When service times are independent of initial 
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lead times, the lead time profile can be expressed in terms of a convolution 
of two measures on R. If initial lead times equal zero, the lead time profile 
provides information about time in queue and sojourn times; see Sections 
3.1 and 3.2. The case in which initial lead times depend linearly on service 
times leads to a heuristic rule for setting "realistic" deadlines for jobs in a 
processor sharing queue; see Section 3.3. 

3.1. Independence of service times and initial lead times. Suppose that 
for each r £lZ, the service times and initial lead times of arriving jobs are 
independent, that is, 

(3.1) e = u r x\ r , 

where v r is the service time distribution on M + and A 7 * is the initial lead 
time distribution on R. Letting A r (C) = X r (rC) for Borel sets CcK, we 
have 



So, by (2.25) 



d r = v r x A r . 



v x A — > v as r — ► oo 



and i) = !/xA, where v and A are the weak limits of v r and A r . By Theo- 
rem 2.2 and the continuous mapping theorem, 

(3.2) Z r (-) ovr" 1 =>.0f*(O ovr" 1 asr^oo. 

In the present setting, t?g = and for z > 0, x E R+ and j/GK, 

(3-3) (l(x,oo)x(y,oo)>^e) = a / ^ ((^ + , oo)) A((y + n, oo)) du. 

•/ 

To compute the limiting lead time profile, set x = in (3.3) and observe 
that #*({0} X R) = 0. Then for y £ R and z > 0, 



(1r + x(j/,oo),0 = « / KC" 2 1 ,oo))A((?/ + n,oo))d'u. 

*/ 

Thus, using Definition 2.1, 

(lK + x(-oo,y])^e) = (^e) ~ x (j/,oo) j $e) 

(3.4) 

/■oo 

= z — a v{(uz~ l , oo))\((y + u, oo)) du. 



Since, by (2.21), 
(3.5) z = za v ((s, oo)) ds = a / ^((mz" 1 , oo)) du, 
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we can deduce from (3.4) that 



oo 



(3-6) (lR + x(-oo,j/],^e) =d I v((uz , oo))A((-oo, y + it]) du 

J 

Letting i>* be the Borel measure on R with density 






Note that if v e denotes the excess lifetime distribution of u, then v e has 
density f*{—uz) for That is, P| is related to u e by a scaling factor z 

and a reflection about zero. 

Letting be the zero measure on R, (3.9) holds for all z > 0. Combining 
(3.2) and (3.9) yields an approximation result for the diffusion scaled lead 
time profiles. As r — > oo, 



As a special case of (3.1), suppose, in addition, that the service time 
distribution v r is exponential for each r G 1Z. Then v is exponential with 
mean a" 1 . For z > 0, (3.7) takes the form 



So, by (3.10), when Z*(t) > 0, the limiting lead time profile A * u e at 
time t is equal to the convolution of A with an exponential measure with 
parameter a/Z*(t) and total mass Z*(t) that is reflected about zero. This 
result applies, in particular, to a sequence of M/M/l queues with lead times 
satisfying A r = A for all r. Thus, it verifies a conjecture of Lehoczky [18], 
who supported this statement with an informal analysis and Monte Carlo 
simulations. 

3.2. Time in queue, sojourn times and the snapshot. An alternative spe- 
cial case of (3.1) can be used to study the profile of times in queue for the 
processor sharing discipline. A job's time in queue at time t is the amount 
of time spent in the buffer by time t. On diffusion scale, the profile of times 
in queue f r (t) at time t is a random finite Borel measure on R + . It can be 



(3.10) 





u<0 
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directly deduced from the diffusion scaled lead time profile by assuming that 
initial lead times are zero for all r £ 1Z. A job's time in queue at time t then 
equals the absolute value of its current lead time; for y > 0, 

(3- 11 ) (l[j/,oo)>^ r (-)) = (1r+x(-oo ,-y], Zr (-))- 

Since (3.1) holds when A r = So for all r 6 7Z, the limiting lead time profile 

■de tt^ 1 is given by (3.9) with A = 5o; for z = 0, o ir~ l equals the zero 
measure on E and for z > and y € R, a short computation using (3.9) and 
(3.7) yields 

(3.12) (1r+x (-00,^,0 = {zv e ([- yz -\ oo)), y< ! 

Thus, if t*(-) denotes the limiting profile of times in queue, then for t > 
and y > 0, 

^.idj UM- T W/-| 0) Z*(t) = 0. 

A related computation reveals the limiting sojourn time behavior, which 
is informally described below. A job's sojourn time is the time between 
its arrival at the buffer and its service completion. Consider the empirical 
distribution of a collection of sojourn times, sampled over a time interval 
that is small on diffusion scale. This distribution can be studied as the 
interval varies over diffusion scaled time. Since the interval will be taken to 
be asymptotically small, it can be defined in terms of residual service times. 

Fix T > 1 and t £ [0, T) such that Z*(t) > 0. For e > and reK, let J r £ (t) 
denote the set of jobs with residual service time in (0, e] at (diffusion scaled) 
time t. Let fj(t) be the diffusion scaled distribution of times in queue for 
those jobs in JJ(t). Once again, assume that A r = 5$ for all r so that a job's 
time in queue is the absolute value of its lead time. Then for y > 0, 

(1[0,i/)>t£(*)) = ¥?((l[o, e ]xR> '( i )))(l[0,e]x(-y,oo), & '(*)) > 

where ip(x) = 1/x for x > 0, with <p(0) = 0. Let <r£(i) denote the diffusion 
scaled distribution of sojourn times for those jobs in JJ(t). That is, if cj-* is 
the unsealed sojourn time of job i € Je(t), then for y > 0, 

(3.14) (l {0>y) ,al(t)) = V ((l [0i£]xR ,Z r (t))) J2 W^O- 

For small e, the probability measure Tj(i) approximates ££(4) in the 
following sense: if job i € Jg(t), then in the unsealed system, vl(r 2 t) < e. 
Since the sequence {Z r (-)} is tight, there exists C > such that, with high 
probability, sup sg r T i Z T (s) < C. Assume that r is sufficiently large that 
r 2 t + reC < r 2 T. Then, with high probability, job i receives service at rate 
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at least (rC) _1 if it is present in the buffer during the entire time interval 
[r 2 t,r 2 t + reC] C [0,r 2 T]. It must therefore complete service at or before 
time r 2 t + reC. Since its time in queue increases by at most reC before it 
departs, its sojourn time o~\ will differ from its time in queue at time r 2 t by 
at most reC. So, by (3.14), its diffusion scaled sojourn time will differ from 
its diffusion scaled time in queue at time t by at most eC. Taking e to be 
small compared to C, it is easy to show that rj(i) and (J r e {t) are close in the 
weak topology. As r — > oo, the number of jobs in J£(t) grows without bound 
for all e > 0, yielding nontrivial distributions fj(i) and cr r £ (t). 

By letting e — > after r — > oo, one obtains a limiting distribution cr*(t). 
By Theorem 2.2 and the continuous mapping theorem, rj(i) can be approx- 
imated by T*(t), where, for y>0, 

(l [0>J/ ),T*(t)) = (l [0)£ ]x(- 2/ ,oc) > C (t) )v((l[0 >£ ]x R ^f (t) »- 

For z > and y > 0, a short computation using (3.3) with X = So yields 

(l[o, e ]x(-y,oo),flg) = ^([O.yz" 1 )) - ^ e ([e,g + j/z' 1 )) 

(l[0,e]xR,^e) l/ e ([0,Oo)) - I/ e ([e,Oo)) 

= z^([0,e)) -Uedyz^^ + yz" 1 )) 
f e ([0,e)) 

Thus, for z > and y > such that v({yz~ 1 }) = 0, 

,. (l[0,e]x(-j/,oo))^e) /r -1 u / Fn -Inn 

(1[0, E ]XR>^) 
Conclude that for y > 0, 

C 7*(t)([0,y))=K[0,y^(t)" 1 )). 

Since ^g(t) and fj(t) are close for small e, the measure cr*(t) may be 
interpreted as a limiting sojourn time distribution of jobs "departing the 
system at diffusion scaled time tP Note that the order in which limits are 
taken here is the only way to obtain a nontrivial limit; if one takes e — > in 
the prelimit, then fj(t) converges to the zero measure because Z r (t) does 
not charge {0} x R. Also, note that the limit cr*(t) has been obtained for 
fixed t when Z*(t) > 0, not as a process. 

An implication of (3.14) and the subsequent computations is that sojourn 
times in the rth system are of order r. Compare this to the state descrip- 
tor: since Z r {-) = r~ 1 Z r (r 2 -) converges to a process with continuous sample 
paths, the state descriptor varies slowly on the time scale of sojourn times 
when r is large. That is, jobs observe only very small changes in the system 
state, including the queue length, during their stay in the system. This is 
analogous to Reiman's snapshot principle (see [18, 22]). 
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3.3. Linear dependence of initial lead times and service times. By the 
results in the previous section, the amount of time a job spends in a heavily 
loaded single server processor sharing queue is roughly proportional to its 
service time and the queue length observed upon arrival. Suppose that in 
the rth system, jobs have "realistic" deadlines, that is, l\ = crv\ for all 
i = 1,2, ... , where c> is a constant corresponding to the (diffusion scaled) 
queue length a job hopes to find upon arrival. Then for x £ R+ and 

(![a;,oo)x[i/,t»)^) =V{[X Vj/C -1 ,Oo)). 

By Theorem 2.2, Z r {-) tff * ( ' } , where for z > 0, x £ R + and y £ R, 

poo 

(l[x,oo)x[j/,oo)) $e) = za / (1[(z+u)V(s/+zm)c~\oo)> zy ) <^ n - 



Setting x = yields, after some computation, 



^ R + x ^°°)' e/ \ci/ e ([yc- 1 ,oo))-(c-z)i/ e ([y(c-z)- 1 ,co)), y > 0, 

for < z < c, where the last term above is understood to be zero when z = c. 
For z > c, 

1 1 ,o*\ _ /2 + (c-z)i/ e ([y(c-2;)- 1 ,cx3)), y < 0, 

(lR + x ^'^-lci/ e (bc- 1 ,oo)) ) y>0. 

As the snapshot principle predicts, there is no lateness in the limiting system 
when Z*(t) < c and therefore almost no lateness in the prelimit when r is 
large and Z r (t) <rc + o(r). If Z*(t) > c [r is large and Z r (t) > r(c + e) for 
some e > 0], then jobs become late before they exit the system. Thus, c acts 
as a threshold for Z*{-), governing when the system exhibits lateness. Since 
the law of Z*(-) is known [9], one can influence the proportion of missed 
deadlines by choosing c appropriately. For example, in cases where 7 > 
(corresponding to a sequence of subcritical prelimit models), Z*(-) has a 
steady state distribution [9]; by choosing c as a particular quantile of this 
distribution, one could achieve a desired long-run average lateness rate. 

4. Tightness. The remainder of the paper is devoted to proving Theorem 
2.2. The goal of this section is to prove that 7£-indexed sequences of sample 
paths chosen from {Z r,m (-):m < [rT\} are relatively compact with high 
probability. 

4.1. Preliminary results. Assumption (A) has some important conse- 
quences which are reviewed below for later reference. Recall that W r (-) = 
{x,Z r (-)) is the diffusion scaled workload process for the rth model. Let 
W*{-) be a reflected Brownian motion on R + with drift —7, variance a(a 2 + 
b 2 ) and initial value iy*(0) equal in distribution to G). 
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Proposition 4.1. As r — > oo, the sequence of diffusion scaled workload 
processes {(x,Z r (-))} converges in distribution to W*{-). 

Proof. It is well known that (A) is sufficient to imply the above result 
for the workload process of any single-server, single-class queue operating 
under a work conserving service discipline (including the present processor 
sharing model). The use of functional central limit theorems and continuous 
mappings to prove such results goes back to [11]. See [24] for a detailed 
account covering a more general setting than the present one. □ 

A similar result holds for the diffusion scaled queue length process Z r (-) = 
(l,Z r (-)). Assume (A) and let W*(-) be the reflected Brownian motion spec- 
ified above. Let C$ = 2a(l + aV) -1 and define Z*(-) = C#W*{-). 

Proposition 4.2. As r — > oo, the sequence of diffusion scaled queue 
length processes {Z r (-)} converges in distribution to Z*{-). 

PROOF. See Theorem 2.3 in [9] and Corollary 2.4 in [9]. □ 

A simple corollary will also be needed. Assume (A) and let T > 1 and 
£,r,e(0,l). 

Corollary 4.3. For all L>1, 

(4.1) liminfP'Y sup \(l,Z r ' m {t)} - (l,Z r ' m (0))| < e) > 1 - rj. 

te[o,L] 

Furthermore, there exists M > 1 such that for all L > 1, 

(4.2) liminfP r f sup (l,Z r > m (t)} V ( X ,Z r ' m (t)) < m] > 1 - 77. 

JWO ° \m<[rT\, J 

te[o,L] 

Proof. Expand the definition of Z r ' m (-) to rewrite in terms of diffusion 
scaling. It suffices to show that for all L > 1, 

liminfP r f sup \Z r {t + h)- Z r {t)\ < e] > 1 - 77, 
r ^°° V te[o,T], J 

he^Lr- 1 ] 

and that there exists M > 1 such that for all L > 1, 

liminfP'Y sup Z r {t) \/ ( X ,Z r (t)) <m) > 1 - 77. 
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Both statements follow from Propositions 4.1 and 4.2; the sequences {{x, 2 r (-))} 
and {Z r (-)} are tight and the limiting processes W*(-) and Z*(-) are a.s. 
continuous. Note that M can indeed be chosen independent of L. □ 

The following result concerning the invariant manifold will be needed. 

Definition 4.4 (Lifting map). For a probability measure i? on H + sat- 
isfying (2.19)-(2.21), let A# :M. + — > M,? be the associated lifting map given 
by 

(4.3) A*z = 0S, z£R + . 

Lemma 4.5. XTie map A^:M + — > M$ zs continuous and M# C M is 
dosed. 

PROOF. By (2.21), (2.29) and a change of variables, = y for all 

y > 0. Let {z n } C M+ satisfy z n -> z as n -> oo. Then (l,^ n ) -» <1, If 
z = 0, then t?g n — > 0. Assume without loss of generality that z, z n > for 
all n and let B C H+ be closed. By (2.29) and two applications of Fatou's 
lemma, 

limsup^g n (B) <a ^limsupl B+ ( uz -i iU yv) du. 
Since B is closed, 

POO 

\\msu V r™{B)<a / {l B+{uz -i )U) ,$) du = #* e (B). 

n— >oo JO 

Conclude from the Portmanteau theorem that z" 1 ^™ — > z~ lr &l and there- 
fore, For the second statement, let {$e n } C M# satisfy -^-> £ 
as n-> oo. Then z n — ► Since A$ is continuous, £ = i?e £ M$. □ 

4.2. Functional Glivenko-Cantelli estimate. Proofs appearing below will 
frequently rely on a functional uniform Glivenko-Cantelli estimate for the 
measure valued arrival process specified in Section 2.1. They require uni- 
form control over the empirical distributions of random vectors (i>[,Z[r -1 ) 
associated with jobs arriving during an interval [s,t]. More explicitly, for 
all intervals [s,t] with length bounded by a fixed constant L, the random 
measure 

rE r ' m {t) 

pr,m -I \ <- 

L s,t -~ 1^ d H,lrr- 1 ) 
i=rE r > m (s)+l 
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must be asymptotically approximated, with high probability, by the measure 
a r (t — s)t? r . For fixed m, [s,t] and for a fixed Borel set B C H+, it follows 
from (A) and the weak law of large numbers that 

(4.4) £ r s f(B)-a r (t-s)$ r (B)^0 asr^oo. 

The present setting requires this approximation to be uniform over a suffi- 
ciently rich class of subsets of H_|_ and, for each T,L> 1, uniform over 
all m < [rT\ and all intervals [s,t] C [0,L]. This section addresses this issue 
using well-known results from empirical process theory. 

Let f be a family of Borel-measurable functions / : H + — > R and let 
F :M + — > R be a Borel-measurable envelope for Y, that is, |/| < F for all 
/ G y. For f £ y and a Borel probability measure Q on H+, let ||/||q )2 = 
(/ f 2 dQ) 1/2 . Let l°°(Y) be the space of bounded functions Q : Y -> R, equipped 
with the supremum norm = sup^ e y |<5(/)|. For a random measure 

£6M (defined on fT) satisfying (F, £) < oo almost surely, identify £ with 
a random element C £ i 00 ^) satisfying £(/) = (/, £) almost surely. This 
identification is unique in law. Then for r £ TZ, m < [rT\ and [s,t] C [0,L], 
the random measure C r s '" 1 — a r (t — s)$ r also represents a random element of 
£°°(Y). A uniform version of (4.4) requires that for an appropriate family 

Y of indicator functions, Clf^ — a r (t — s)3 r — > in £°°('^ / ), uniformly in 
m< \ tT\ and [s,t] C [0,L]. 

This uniform convergence holds under easily verifiable conditions on Y . 
The first main condition guarantees a minimal amount of measurability 
when dealing with suprema over "V . Call a family Y a Borel-measurable 
class if for each n £ N and (ei, . . . , e n ) £ { — 1, l} n , the map 



7-' 

is Borel measurable on H™ . The condition states that for all 5 > and r £lZ, 
the families ^ = {/ - 5 : f,g ef, \\f-g\\^ 2 < 5} and = {(/ - 5 ) 2 : /, 9 G 
^} are Borel-measurable classes. 

The second main condition ensures that Y is not too large, as mea- 
sured by the following combinatorial property. Given a set S, a Vapnik- 
Cervonenkis class ( VC-class) is a collection 5? of subsets of S with finite 
Vapnik-Cervonenkis index ( VC-index) . The VC-index of a collection of sub- 
sets is the smallest integer n (oo if none exists) such that 5? shatters no re- 
point subset {x\, . . . , x n } C S. A collection ^ shatters {x±, . . . , :r n } if each of 
the 2 n subsets of {x\, . . . , x n } can be picked out by J? 7 , that is, can be written 
A Pi {x\ , . . . , x n } for some A £ 5? . Finally, a family of functions ^ : S — > R 
is a VC-class if the collection of subgraphs {{(x,y) :y < f(x)} : / £ ^} is a 
VC-class of sets in 5 x R. For a more detailed discussion of VC-classes, see, 
for example, [23], Chapter 2.6. 
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VC-classes satisfy the uniform entropy condition: Let Y be a VC-class 
with envelope F. Let =2 be the set of finite discrete probability measures Q 
on H + such that ||.F||q2 > 0. For Q G J and a Borel- measurable function 
f:M + ^R satisfying ||>|| Q)2 < oo, let B f {e) = {g G V : \\j - g\\ Q , 2 < e} de- 
note the L2((5)-ball in "V with center / and radius e. Let N(e, "V ,L 2 (Q)) be 
the smallest number of balls Bf{e) needed to cover Y . Then ~f satisfies 

(4.5) [°° sup JlogN(e\\F\\ Qt2 ,y, L 2 (Q)) de < oo. 

Jo Qe£ v 

See Definition 2.1.5, (2.5.1) and Theorem 2.6.7 in [23]. 

The uniform entropy condition (4.5) is used in the next lemma to show 
that, for large k G N, the £°°('^ / )-norm of the random measure 

M) ^ = 7^E(%^-i)-^) 

satisfies a tail bound that is uniform in r G 1Z. In particular, (4.5) implies 
that for each r ElZ, the class "V admits a & ' -Brownian bridge, that is, a tight, 
Borel- measurable version of the zero-mean Gaussian process {Gl (/) : / G ^} 
with covariance function 

v r [SlU)S r M] = (fgJl - (fJ r )(gJ r )- 

Moreover, the expected norms E r [|| (/, Ql) \\ y] are uniformly bounded in r 
and for large k, the random measures {Q^'-r G 1Z} are uniformly close in 
to {Gl'.r G 1Z}. This quickly leads to the desired tail bound; see 
Lemma 4.6 below. Finally, Lemma 4.7 below employs this tail bound to 
establish the functional Glivenko-Cantelli estimate for Cf^ '. 

The proofs in this section use outer probabilities and outer expectations 
to avoid cumbersome verifications of the measurability of certain suprema 
over uncountable sets. For each r G 1Z and each subset Scfl r , define the 
outer probability 

P r (S) = M{P r (B) :ScB,BeF} 

and the outer expectation of a map X : — > M by 

E r [X] = inf{E r [Y] : X < Y, Y is ^-measurable}. 

Assume (A) and let "V be a VC-class of Borel- measurable functions / : H+ — 
]R with Borel-measurable envelope F such that and are Borel- 
measurable classes for all r G 1Z and 5 > and such that 

(4.7) lim sup(F 2 l {F>N} J r )=0. 

N—foa r g7^ 



24 H. C. GROMOLL AND L. KRUK 

Lemma 4.6. For all q > 1 and x > 2, there exist M < oo and ko £ N 
swc/i £/iai k> ko implies 

supP r (\\(f,Cl)\\ r >x)<^. 

The constant M does not depend on x. 

Proof. Fix q > 1 and x > 2. For Q E define 

(4.8) M0) = (IW)lk-z + l) + Al. 
Then 

(4.9) supP'(||(/,^>|| r > x) < supE^M^)]. 

The function /i x is an element of BLi(^), the set of all h:H,°°{Y) — > R sat- 
isfying [|/i||oo < 1 and |/i(e?i) - < - g 2 \\f for all £ £°°(f). 
Since "V is a VC-class, it satisfies the uniform entropy condition (4.5). More- 
over, the envelope F satisfies (4.7) and "V^,^ satisfy the indicated mea- 
surability conditions. Thus, "V satisfies the assumptions of Theorem 2.8.3 in 
[23] which asserts that "V is pre- Gaussian and Donsker, uniformly in r £ 7Z 
(see Section 2.8.2 in [23]). In particular, ~f satisfies 

(4.10) supEl|(/,Sr)|H<oo 
and since h x £ BLi, 

(4.11) lim sup lE'^C^)] - E r [h x (gi)}\ = 0. 

By (4.9) and (4.11), there exists ko £ N such that k > ko implies 

S npt> r (\\(f,g r k )\\ r >x)< S u V E r [h x (g;)] + x -«. 

Apply (4.8) and Markov's inequality to the right-hand side to obtain 
su P P r (|K/,^)|| r > x) < supP'XIiaSDIk >x- 1) + 2T« 

<x-^2^upB r [\\(f,gi)\\^] + l 

Let M be the last term in parentheses, which does not depend on x. For each 
r € 71, the Brownian bridge is separable and Gaussian with ||(/, 
finite almost surely. Thus, there exists a constant C such that for all r £ 7Z, 

E r [\\(f,g:)\\^]<CE r [\\(f,gi)\\^ 

(see Proposition A. 2. 4 in [23]). Conclude from (4.10) that M < oo. □ 
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The next lemma establishes the functional Glivenko-Cantelli estimate for 
£ r s '™. Assume (A) and let T,L > 1 and e G (0, 1). Let V be a VC-class of 
Borel- measurable functions / : H_|_ — > R such that and "V^ are Borel- 
measurable classes for all r G 7Z and 5 > and such that (4.7) holds. 

Lemma 4.7. There exist events 

njc/sup sup |(/,^;n-a r (<-«)(/,^)|<e}. 

l/Gr m<|rTJ, J 
[«,t]C[0,L] 

smc/i i/taf limf-^oo P r (QJ) = 1. 



Proof. For each r€TZ, let 

(4.12) V = \ sup ||(/,^;n-a r (*-a)(/,^>|k >e}- 

I m<\rT\, ) 
[s,t]c[0,L] 

It suffices to show that 

(4.13) lim P r (T r ) = 0. 
Since 

{f,C:f)-a r (t-s){fJ r ) = ((f,£ r (m + t))-a r (m + t){f,$ r )) 

-({f,C-(m + s))-a r (m + s){f,3 r )), 
the supremum in (4.12) is bounded above by 

sup \\((f,C r (t))-a r t(f,^})-((f,C r (s)}-a r s(f,^})\\ r . 

[s,t]c[0,rT+L], 
\t-s\<L 

Rewrite this as 

sup \\af,C r (rt))-a r rt(f,$ r )) - ((f,C r (rs))-a r rs(f,^))\\y 

[s,t]c[0,T+L/r], 
\t-s\<L/r 

For u > 0, let ¥ \u) = r- 2 E r (r 2 u) and E r (u) = E r (ru) - a r ru. Then 

rE r (ru) 

(f,£ r (ru))-a r ru(f,fr) = - £ (/(<> ^) ~ </>^» 



r 

i=i 



+ (£ r (ru)-c/ru)(/,?? r ) 
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Introducing the notation X T k {f) = {k/r 2 ) l l 2 (f,G r k ) for k E N and / e f , the 
supremum in (4.12) is bounded above by 

sup \\JT -, ( Jf)-Xy r Jf) + (E r (t)-E r (s))(fJ r )\W- 

[s,t]c[0,T+L/r], W { > 

\t—s\<L/r 

Let 5 > 0. Deduce from the previous remarks and the inequalities || (/, $ r ) \\y < 
(F, ■# r ) and limsnp r ^ 00 L / r < T A 5 that 



lim P(T r ) <limsu P P r sup \E r (t) - E r (s)\(F,$ r ) > - 

r^oo r^oo \s,t<2T, 2 

\t-s\<8 

(4.14) 

+ iimsu P p r ( sup \\xy E r (t) (f)-xy B r {s) (f)\\ y >^ 



s,t<2T, 
t~s\<8 



By (2.22)-(2.24) and the functional central limit theorem for renewal pro- 
cesses, the sequence {E r (-)} converges in distribution to a Brownian motion 
which is almost surely uniformly continuous on compact time intervals [2]. 
Since sup re ^(-F, ■& T ) < oo by (4.7), we can conclude that the first right-hand 
term in (4.14) converges to zero as 5 — > 0. To show the same for the second 
term, the proof of Theorem 2.12.1 in [23] has been adapted. 
For a, b £ [0, oo) , define =Nn [a, b] . It suffices to show that 



(4.15) limsupPM max 



S u P \\x; 2¥{t) (f)-xy r{s) (f)\\ r >- 



1 s,te[(j-l)8,jS] 

converges to zero as 5 — > 0. For each r E 1Z, bound the outer probability of 
the maximum by the sum of outer probabilities and discretize the time index 
in the supremum. Then (4.15) is bounded above by 

[2T/6] 



(4.16) limsup J2 



max \\X[(f)-X r k (f)\\ r >- 

r 2 £ r ((j-l),5) 



Let fig = {sup ug [ 0) 2T+<5] 1-^ ( u ) ~ ot r u\ > a5}. Then on the complement of S7q , 
for each j = 1, . . . , \2T/5\ , 

(a r jS -a5) + < W(j8) < a r j5 + at6. 
Deduce that (4.16) is bounded above by 

\2T/5] 

lim sup P r 

r— >oo , 



max 

r A (a r {j — 1)5 — aS)^ 



r r 



£ 



+ lim sup 



"2T" 
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The second term equals zero by (2.22)-(2.24) and the functional weak law 
of large numbers for renewal processes. Since Q r k has stationary increments 
as a function of k and limsup r ^ 00 Q r < 2a, the previous term is bounded 
above by 

Vk 



(4.17) 



lim sup 



2T 



max 

k<[r 2 4a5\ 



£ 

^4 



By Ottaviani's inequality (see Proposition A. 1.1 in [23]) and the stationary 
increments of Q r . , (4.17) is bounded above by 



(4.18) 



Assume that 5 is sufficiently small that e{lQya5) 1 > 2. By Lemma 4.6, 
there exist M and ko £ N such that k > ko implies 

\\{f,Gk)h> — = i < 



[2T/51 (||(/, g[ r24Q(5J ) || r > e/(16 Va*)) 

lim sup : - =— . 

™ 1 - max fe < Lr24a5J P r (\\(f,G r k )\\y > er/(8y/k)) 



(4.19) 



sup P T 



16v a5 J 



M. 



Since [^ 2 4a5J — > oo as r — > oo, the limit superior of the numerator in (4.18) 
is bounded above by |~2T/<5] (16\/a5/e) 3 M, which can be made arbitrarily 
small by choosing 5 sufficiently close to 0. It remains to show that the limit 
inferior, as r — ► oo, of the denominator in (4.18) is bounded away from zero 
as 5 -► 0. Observe that er^^/Pr 1 > eiwV^S)- 1 for k < [r 2 Aa5\ . So, by 
(4.19), the terms in the maximum indexed by k > ko are bounded above 
(uniformly in r) by (16\/a5 / e) 3 M , which can be made arbitrarily small by 
choosing 5 sufficiently close to 0. For k < ko, 



P r \\{f,Gi)\W> 



er 
8Vk 



<P r \\(f,Q r k )h> 



er 
8JE 



which converges to zero as r — > oo. Conclude that (4.18) converges to zero 
as 5^0. □ 



It remains to define a class of subsets of H + that is sufficiently rich to 
characterize elements of M, yet sufficiently small to be a VC-class. 



Definition 4.8. Let (a, b) denote the closure in R of an interval (a, b) C 
R and define 

si = {[x, oo) x (y, oo) : x £ [0, oo), y G [— oo, oo)}. 

Note that since the above definition allows y = —oo, the collection ^ 
includes not only subsets of M + that are translations of the right upper 
quadrant, but also all right half-spaces. 
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Lemma 4.9. The family "V = {\a '■ A £ is a VC-class of Borel-measurable 
functions with Borel-measurable envelope F = 1. Moreover, and Y^' are 
Borel-measurable classes for all r £lZ and 5 > and the envelope F satis- 
fies (4.7). 

Proof. Since "V contains only indicator functions, it suffices to show 
that si is a VC-class. Let {a,b,c} C H + be a three-point subset. It is im- 
possible for g/ to pick out all three two-point subsets of {a,b,c}. Thus, g/ 
shatters no three-point subset of H + and so has VC-index at most 3. The 
remaining assertions are evident. □ 

4.3. Dynamic equation. This section introduces the dynamic equation 
satisfied by the measure valued state descriptor. For a subset B C H + and 
w G HL|_, define l~g(w) = (1b, S^) (recall the definition of <5+ from Section 
2.1). For each r G 1Z, the state descriptor of the rth model satisfies the 
following equation almost surely: for each Borel set B C H_|_ and all t,h>0, 

Z r (t + h)(B) = Z r (t)(B + (Sl t+h ,h)) 

E r (t+h) 

+ E i B (vut+h),i:(t+h)). 

i=E r (t)+l 

This follows from (2.3)-(2.8) after some simplification. Applying the fluid 
scaling (2.10)-(2.15) and (2.17), the above equation becomes 

Z r (t + h)(B) = Z r (t)(B + (Sl t+h ,h)) 

(4.20) 

rE r (t+h) 

+ - E i + B m+h)ji(t+h)). 

i=rE r (t)+l 

After rescaling and shifting, the equation takes the following form. For each 
r &1Z, m< \ rT\ and t,h>0, almost surely for each Borel set B C H+, 

Z r > m (t + h)(B) = Z r ' m (t)(B + (Sf£ fcJ h)) 

(4.21) 

rE r ' m {t+h) 

+ - E l B (vr(t + h)X' m (t + h))- 

i=rE r > m (t)+l 

Equations (4.20) and (4.21) are called the dynamic equations for Z r (-) and 
J r ' m (-). 

Subsequent proofs use estimates obtained from (4.21). Two estimates re- 
sult from bounding the summands by 1 and optionally bounding the first 
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term on the right-hand side by its total mass; for each B C H+ and t,h>0, 
Z r > m (t + h)(B) < Z r > m (t)(B + (S r t f +h , h)) + C r t f +h {H + ) 

(4.22) 

<z^(t)(m + ) + £ r t f +h (m + ). 

Two more estimates follow from (4.21) by simply ignoring any arrivals; for 
each B CM + and t,h>0, 

(4.23) Z r > m (t)(B + {S r t f +h , h)) < Z r > m {t + h){B) < Z r > m {t + h)(M + ). 

4.4. Compact containment. The following lemma establishes the com- 
pact containment of the fluid scaled state descriptor on [0, [rT\ +L\. Assume 
(A) and let T > 1 and rj G (0, 1). 

Lemma 4.10. There exists a compact set K C M such that for all L > 1, 
liminf P r (Z r (t) G K for all t G [0, \rT\ + L\) > 1 - n. 

r — >oo 

Proof. Since 7r + :(x,y) i— ► \y\ is continuous, (2.25) and the Skorohod 
representation theorem imply the existence of M+-valued random variables 
X r ~ -d r o tt^ 1 and X ~ $ o -zr^ 1 such that X r — > X almost surely. Conse- 
quently, there exists an M+-valued random variable Y such that 

(4.24) Y = supX r a.s. 

r 

Let /x be the law of Y on R+. Since L 2 (fi) contains unbounded functions, 
there exists a continuous increasing unbounded tp : R + — > M + such that (ip 2 ,fi) < 
oo. This implies that 

<ty o 7T+) V) = B[^X) 2 ] < n^(Yf] < oo. 

Thus, the one-element set "¥ = {ip o ir + } is a VC-class with envelope F = 
tp o 7r + satisfying (4.7). 

Since M is a Polish space, (2.32) and Prohorov's theorem imply that the 
sequence {Z r (Q) : r G 1Z} is tight: there exists a compact Ko C M such that 

(4.25) liminf P r (Z r (0) G K ) > 1 - |. 

By (4.25), Corollary 4.3 and Lemma 4.7, there exist M > 1 and events 
n;<=( sup |(^ -cfty o7r+,^)l <i), ren, 

lm<\_2rT\ ) 
such that for all L > 1, the events 
W 2 = {Z r (0)£K }, 

n r 3 = S sup (l,Z r (t)}V( X ,Z r (t)}<M\, 

Ue[0,LrTJ+L] J 

Q5 = n Q r 2 n n§ 
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satisfy 

(4.26) liminfP r (^) > 1 - n. 

r — >oo 

Assume henceforth in the proof that all random objects are evaluated at 
some outcome wsfij. 

Since Ko is compact, there exists a sequence of real numbers {x k : k S N} 
such that for all k, 

(4.27) sup e(([0,s fc ] x [-x k ,x k ]) c ) < -. 

For each k G N, choose > such that i^(y k ) > A; 2 and define 

4 = [o,fc], 

</fc = [-{yk + k),y k }, 
Rk = h x Jfc. 
Consider a sequence — > as k — > oo. The set 

(4.28) K = {£ G M : (1, £) < M and £(#£) < a fc for all A;} 

is precompact ([12], Theorem A 7.5). Thus, by (4.26) and the definition of 
Jig, it suffices to identify a null sequence {a^} such that for each L > 1 and 
each r > L, 

sup Z r {t){R c k )<a k for all A;. 

t€[0,|rT|+L] 

Fix L > 1, t G [0, [rT\ + L] and fe G N. Assume that r > L. Observe that 

(4.29) Z r (t)(R c k )=Z r (t)(I c k x R) + Z r (t)(I k x J c k ). 

Recall that x( x iU) =x. By Markov's inequality and the definition of fijj, 

(4.30) Z r (t){It x R) = Z r (t)((k, oo) x R) < ±<x,Z r (t)> < 

Let AA = {sG [(*- fe)+,t]:(l,Z r (s)) = 0}; if M^0, fix some t e J\f and if 
= 0, let r = (t - k)+. By (4.20), 

z r (t)(4 x r k ) = z r (r)((i k x r k ) + (a^ t ,t - r )) 

(4.31) 

r£T(t) 

+ - E lix^^rwJK*))- 

i=r_E r (T) + l 

If M ^ 0, then 2 r (r) = and 

(4.32) £ r (r)((4x J-) + (5;,,t-r))=0. 
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If J\f = 0, then either t = or r = t - k. Suppose r = 0. Then by the defini- 
tions of y k and J k , by (4.27) and by the inequality t — r < k, 

(4.33) Z r (r)((I k x 4) + (S r Tit ,t- r)) < Z r (0)(R+ x [-y k ,y k ] c ) < ±. 

Suppose r = t — k. Then, by definition of and the fact that M = 0, we 
have 

"* t-r k 



M M 

In this case, Markov's inequality and the definition of yield 
Z r (r)((I k x J c k ) + (S r T>t ,t- r)) < Z r (r)([k/M,oc) x R) 

(4-34) <^ (X) ^(r)> 

- fc ' 

We can then deduce from (4.31)-(4.34) that 

M 2 1 rS? ' (t) 

(4.35) z r (t)(i k x4)< — + - e ij XJ o(«r(<),rr(<))- 

i=r\E r (Y)+l 

For each i appearing in the sum, (2.5) and (2.15) imply that 

Tl(t) = l\r~ x - t + U[ r" 1 > l\r- x -t + T> ^r" 1 - k. 
Thus, ll(t) £ [-{y k + k),y k ] c implies that Z[r _1 G [-yfc,?/fc] c , yielding 

ljx^l(«(*)) < 
By definition of ^ and 

Vxb.rfK.^" 1 ) < ^)^(l^ r_1 l) ^ pV-diTr- 1 !). 

Deduce from (4.35) that 

M 2 1 1 
^(f)(/*xJ£)< — + - E ^(l^- 1 !). 

i=r_B r (r)+l 

Bounding and rewriting the sum, we obtain 

M 2 1 L<J 1 r ^ r ' m(1) 

z r (t)(i k x j c k )< — + ^ E - E VKI^- 1 !). 

m=|rj i=r,E r > m (0)+l 
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By assumption, t < rT + L = r(T + Lr 1 ) < 2rT. Rewriting the sum again 
and using the definition of we get 

M 2 1 LiJ 
Z' r (t)(I k xJi)<— + - £ (^ovr + ,4;n 

m= |_rj 

(4.36) 

. M 2 fe + 1 . r . 1 , 

<-^ + -p-(« (^ 07r +^ ) + !)■ 

By (2.22), (4.24) and the definition of ^, there exists K < oo such that 

(4.37) supc/(V>o 7r+,r) = supa r E[V>pf r )] < sup a r E[V>(^)] < if. 

rG72. r&TZ r£lZ 

Combining (4.29), (4.30), (4.36) and (4.37) yields 

(4.38) Z r (t)(Rl) < \{M + M 2 + 2(K + 1)). 

k 

Let afc equal the right-hand side of (4.38). Since a k does not depend on 
L and since (4.38) holds on % for all L > 1, r > L and t £ [0, [rT] + L], the 
proof is complete. □ 

4.5. Asymptotic regularity. This section establishes that, as r — ► oo the 
fluid scaled state descriptors Z r (-) assign arbitrarily small mass to the 
boundaries of sets A 6 srf ' . This is phrased in terms of ^-enlargements of 
the boundaries of these sets. For B C H + and k > 0, let 9# denote the 
boundary of B and recall that dg = {w £ EI + :mi z& Q B \\w — z\\ < k} is the 
^-enlargement in M + of its boundary. 

The result is established first for the initial condition Z r (0). Assume (A) 
and let e,rj € (0, 1). 

Lemma 4.11. There exists k > such that 

(4.39) liminfP'Y sup J r (0)(^) < e) > 1 - n. 

r->oo y A6jaf y 

Proof. Note that the event in (4.39) is P r -measurable for each r; since 
the random Borel measure Z r (0) is continuous from below, the event can 
be rewritten using the supremum over a countable family {c^ }, where q is 
rational and {A n } C have boundary points with rational coordinates. 

By the first component of (2.32) and the Skorohod representation the- 
orem, there exist random measures A r ~ Z r (0) and A ~ G, defined on a 
common probability space (J7,^",P), such that A r A almost surely. It 
suffices to show (4.39) for A r in place of Z r (0). Since M is Polish, the se- 
quence {A r } is strongly tight: there exists a compact K C M such that 

= Hreili^ G K l satisfies 

(4.40) p(n ) >i-v 
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(see [17], Corollary 2). Let K# = K n M#. By (2.30) and the fact that K is 
compact, 

(4.41) A G K,? a.s. on Q . 

Choose a compact C C H+ such that 

su P e(c c )<^ 

(see [12], Theorem A 7.5). For each k > 0, let F be the collection of all 
vertical strips [in, (i + 1)k] x R, i = 0, 1, . . . , and all horizontal strips x 
[jk, (j + l)/c], j G Z. Define 

I" ={/GF:/nC/0}, 

which is finite because C is compact. Note that for each A £ g/ , the set 
d^nC is contained in the union of at most six strips in Ig, . For all k > 
and £ G K, 

S u P e(5i)< su P e(5inc) + ^ 

<6max£(I) + -. 

Consequently, 

liminfpf sup A r (d K A ) < e] > liminfpf imaxA r (J) < — 1 nfi ( 
Apply Fatou's lemma to the right-hand side to obtain 



liminfP sup A r (d K A ) <e 

(4-42) 

> P( I limsupmaxA r (I) < — I n Q ( 
\ l. r— >oo /eig v ' 12 J 

Each I G I£< is closed, so the Portmanteau theorem implies that 

(4.43) limsupmax A r (I) < max A(I) a.s. 

Combining (4.42) and (4.43) yields 

liminfP^ sup A r (d%) < e^j > pQmaxA(I) < n ft ( 

Thus, by (4.40) and (4.41), it remains to show that 

(4.44) lim sup max£(I)=0. 
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If (4.44) fails, then there exist k„ -> 0, {£„} C K#, {I n : I n G l£ n } and 5 > 0, 
such that 

inf ><5 

n 

and such that are either all vertical or all horizontal. Suppose that they 
are all vertical. For each n, choose x n G I n n {R+ x {0}}. Since C and K# 
are compact, assume (by passing to a subsequence if necessary) that x n — ► x 
and £ n — ► £. For to > 0, let J™ = [(x — w) + , x lnjxi Then for each w > 
and n sufficiently large, I n C I™. Thus, 

Kmiiif£ n (0 ><5. 

n—i-oo 

Since I™ is closed, the Portmanteau theorem implies that for each w > 0, 

^)>lkkf^)>i. 

Letting I x = C\w>o^x = i x } x ^> deduce that 

(4.45) e(4) > 5. 

However, £ G M$ because is closed (Lemma 4.5). So £ = i?* for some 
z > 0; by Definition 2.1, 

(4.46) £(Ix) = a ^({x + uz- 1 } xR)du. 



Since $(• x R) is a probability measure on R + , it has at most countably 
many atoms. Thus, 

(4.47) £(4) = 0, 

contradicting (4.45). The argument is identical if {I n } are all horizontal. 
Conclude that (4.44) must hold. □ 

The regularity result is now shown for the entire state descriptor Z r (-). 
Assume (A) and let T, L > 1 and e, rj G (0, 1). 

Lemma 4.12. There exists k > suc/i f/iai 

(4.48) liminfP r fsup sup Z r {t){d%) < £ J > 1 - rj. 

rwoo \A&efte[0,lrTJ+£] / 

Proof. By Lemmas 4.11 and 4.10, there exist kq > and compact K C 
M such that the events 

n; = {su P z r (o)(d%)< £ - 

n r 2 = {Z r (t) G K for all t G [0, \rT\ + L}} 
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satisfy 

(4.49) liminf P r («£ n ST 2 ) > 1 - |. 

Let Rk = G H + : ||w|| < X} for each if > 0. Since K is compact, there 
exist M G (1, oo) and K < oo such that 

(4.50) sup(l,£)<M, 

(4.51) supe(j R^)<£ 

£eK 2 

(see, e.g., [12], Theorem A 7.5). Define a* = sup rg 7^a r , which is finite by 
(2.22). Let h = e^a*)' 1 A (L - l)/2, let k = k A h{2M)~ l and let 5 = 
e^KMh' 1 ])' 1 A2" 1 ). By (4.49) and Lemmas 4.7 and 4.9, there exist 
events 

fiScjsup sup \£ r s f(A)-a r (t-s)^ r (A)\<s\, reK, 

iAea' m<\rT\, ' J 
[s,t]c[0,L] 

such that the events Qq = Q\ n ^ ^3 satisfy 

(4.52) liminfP r (OI) > 1-77. 

r — >oo 

Let fi£ denote the event in (4.48). By (4.52), it suffices to show that Qq C fi* 
for all r G 7?.. Let w£S]q. For the remainder of the proof, all random objects 
are evaluated at this u. 

Consider any t G [0, [rT\ + L] and Define 

(4.53) n=sup{s<t:(l,Z r (s)) = 0} 

if the supremum exists and define t\ = otherwise. Let r = t\ V (i — KM). 
The first step is to show that 

(4.54) z r (T)(d% + (S^ t ,t-r))< £ -. 
If r = 0, this follows from the definition of $7^ because 

and because k < kq and srf is closed under positive translation. Suppose r = 
Ti > 0. Then t = m + s for some m < [rT\ and s G (0, L]. By (4.53), there 
exists a sequence {s n } C [0, s] such that (l,Z r,m (s n )) = for all n. By (4.22) 
and the definition of Q3, 

^(rXflX + (S r T ,u t~r))= Z r ^{ S )(d% + (S r T>t , t - r)) 

<^ m ( Sn )(M + ) + £^ s (H + ) 

< a r (s - s n ) + 5 



S r T ,t= f {l,Z r {s)y l ds>K. 
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for all n G N. Letting s n ] s, 

Z r (rm + (S^ t ,t-r))<6< £ -. 

Suppose that t = t — KM. Since (1, Z r (s)) > for all s G (r, t] , the definition 
of f^2 and (4.50) imply that 

r* 

It-KM 

So, by definition of fi£ and (4.51), 

^(r)(,91 + (S;„ t - r)) < £ r (r)(i^) < |. 

The preceding three cases prove (4.54). 

Using the dynamic equation (4.20), we obtain 

Z r (t)(d%)=Z r (T)(d% + (S r T j,t-T)) 

(4.55) 

r£ r (t) 

+- E fymw®)- 

i=rE r {r)+l 

Let / denote the second right-hand term in (4.55). By (4.54), it remains to 
show that / < e/2. Let N=\(t— r)/!" 1 ] and for each n G {0, . . . , N - 1}, let 
t n = r + nh and t n = t n+ \ A t. Then, using the inequality Iqk (•, •) < lg« (•, •), 

N-l , rE r (t n ) 

(4-56) ^<E- E 

n=0 i=r£ ;r( tn ) +1 

Consider n G {0, . . . , N — 1} and i such that U[r~ 1 G (i n , t n ]. Observe that 
(4-57) Sfnj < S[jr r -i t < Sl n t . 

By (2.3), (2.5), (2.14) and (2.15), we have 

(4.58) l ai (vl(t)JUt)) = laj-Ks^ ^t-i^r-^K.^" 1 )- 
So, letting 

A~ = A + (5[n >t - k, t - i n - «) n H+, 
^4 n = A + (Sl n j -\- K,t — t n -\- k), 

An = ^4 n \ A„ ) 
it follows from (4.57) and (4.58) that 

(4.59) la^vimm^uMX^ 1 )- 
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Conclude from (4.56) and (4.59) that 

N-l , rE T {t n ) N-l 
n=0 i=rE r (t n )+l ™=0 

Note that t n — t n < h < L for all n, so C,\ n t n can be rewritten C^'J 1 for some 
m < [rT\ and [s,t] C [0, L]. Since A~,A+ E &/ for all n < N, we can deduce 
from the definition of £l r 3 that 

N-l 

71=0 

By the definitions of a* and ./V and the fact that t — r < KM, 

N-l 

I<a*hJ2 $ r ( A n) + \KMh-^25. 

n=0 

This implies, by choice of 5, that 

N-l 

(4.60) J<a*/i^^(^) + 7- 

n=0 

If ne{0,...,JV-3}, then 

St +1 t +2 > ^Af" 1 > 2^ 

because < (l,J? r (s)) < M for all s 6 (r,t] and because /i > k2M, by defi- 
nition of k. Thus, for all n £ {0, . . . , A r — 3}, 

Deduce that A~ C -A^ +2 for all n G {0, . . . , iV — 3} and, consequently, that 
A n r\A n+ 2 = 0. Thus, since & r is a probability measure, ^1=0 1 ^ r (^2n) 
and En=o 2 ^ 2 ^ ^ r ( j 42n+i) are both bounded by one. Conclude from (4.60) 
that 

K2a*h + ^-, 
4 

which implies, by choice of h, that I <e/2. □ 



4.6. Oscillation bound. This section contains the main oscillation bound 
that constitutes the first of two main ingredients needed to prove tightness 
of the state descriptors. 
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Definition 4.13. For each L > 1, £(•) G D([0, oo), M) and 5 > 0, define 
the modulus of continuity of £(•) on [0,L] as 

(4.61) W£(C(-),5)= sup sup d[((t + h),((t)}. 

te[o,L-S\ he[o,s] 

Denote the modulus of continuity of ("(•) on [0,L) by 

(4.62) wl_(C(-),*)= sup sup d[C(t + /i),C(t)]. 

te[o,L-5) /ie[o,<5] 

Assume (A) and let T,L > 1 and e, 77 e (0, 1). 
Lemma 4.14. T/iere exists 5 > stjc/i that 



(4.63) liminfPM max w L (£ r ' m (-), 6) < e > 1 - 77. 

Proof. By Lemmas 4.12 and 4.7, there exists k G (0, 1) such that for 
each fixed 5 > 0, the events 

n; = ( max sup £ r ' m (i)([0, K ] x R) < f V 

Lm<LrTJ tG [ 0jL ] 4 J 



^2 = 1 max sup £^(11+) < 2aS \, 



i<\rT\ te[0,L-S] 

satisfy 

(4.64) liminfP r (^) > 1-77. 

r — ►oo 

Fix 5 = Ke 2 (8(a V l))" 1 and let Ql be the event in (4.63). By (4.64), it 
suffices to show that J)q C fi* for each r. Fix r £lZ and w G f^. For the 
remainder of the proof, all random objects are evaluated at this uj. Fix 
m < L^TJ , t 6 [0, L - <S], /i 6 [0, <5] and let Bci+ be closed. By definition of 
the Prohorov metric d[-, •], it suffices to show the two inequalities 

(4.65) Z r ' m (t)(B) < Z r ' m (t + h)(B £ ) + e, 

(4.66) Z r ' m (t + h){B) < Z r ' m (t)(B £ ) + e. 
To show (4.65), we use the definition of fi£ to write 

Z r ' m (t)(B)<Z r ' m (t)([0,K] xR) + Z r ' m (t)(Bn((K,oo) x 

(4.67) 

<i + z r > m (t){Bn({ K ,oc) x 
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Let J = {s G [t,t + h]:(l,Z r ' m {s)) < e/4}. Suppose that 1 = 0. Then 
(l,Z r ' m (s)) >e/4 for all se[t,t + h], which implies that 

r t+s - , 47) 

(4.68) \\(S r t ;™ h ,h)\\< / (l,Z r ' m (s))- 1 ds + 5< — + 5<eA K . 

Consequently, w G B n ((k, oo) x R) implies that u; — /i) G -B e and so 

(4.69) 5 n ((«, co)xK)cB £ + (S r t ^ h , h). 
Deduce from (4.67) that 

z r > m (t)(B) <£ + ^ m (i)(i? £ + (3^, fc)). 

Applying (4.23), we obtain 

(4.70) Z r ' m (t)(£) < | + Z r ' rn (t + /i)(£ e ). 

Suppose that 1^0 and let_ r = inf J. Then (1, J r ' m (r)) < e/4, by right 
continuity. If r > t, then (l,Z r ' m (s)) > e/4 for all s G [i,r), so 

(4.71) S[f = r<l,Z r ' ro (fl))- 1 da < 4(T-t) < — < «. 

7t e e 

By (4.67) and (4.71), 

^ r ' m (t)(S) <| + ^ r ' m (t)(( K ,oo) xR) 

<£ + ^(t)((^f,oo)xR). 

Applying (4.23), we obtain 

(4.72) Z r > m (t)(B) < £ + Z^(r)(H + ) < |. 

Therefore, (4.65) follows because either (4.70) or (4.72) holds. 

To show (4.66), we use (4.22) and the definitions of W 2 and 8 to obtain 

Z^it + h)(B) < Z^{t)(B + (Sl : ™ h , h)) + cr t f +h {n + ) 

(4.73) 

<Z^{t){B + {S r t f +h M + \- 

If 1 = 0, then (4.68) implies that B + (S^. h ,h) C B £ . Therefore, (4.73) 
yields 

Z r ' m (t + h){B) < Z r ' m (t){B £ ) + £ 
If I 7^ 0, then by (4.22) and the definitions of an d <5j 

Z^(t + < Z^(r)(H + ) + ^ h (E+) < \ + 2a5 < £ 
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In both cases, (4.66) holds. 

Conclude from (4.65) and (4.66) that 

d[Z r ' m (t),Z r ' m {t + h)]<e. 

Since m < [rT\ , t G [0, L — S] and h G [0, 5] were arbitrary, 

max w L (Z r ' m (-),5) <£, 

m< [rT] 

which implies that □ 

4.7. Precompactness. This section introduces a sequence of events on 
which the shifted fluid scaled state descriptors {Z r,m (-):m < L r ^"J} have 
a desired list of properties. It will be shown that the tail of this sequence 
has arbitrarily high probability and that certain sequences of sample paths 
chosen from the events are precompact. This will be used in the next section 
to construct fluid approximations to the shifted fluid scaled state descriptors. 

Assume (A) and let T > 1. Fix a positive constant M and a compact 
set KcM. For each L > 1, let C L = {{Kj}f =l , {5 k }k%i) be a collection, 
depending on L, where Kj,5k are positive constants with 5 k — > as k — > oo. 

Definition 4.15. For each L > 1 and each r G 1Z and n G N, define 
events 

O r,1 = ( max sup (1 Vx,^ r,m W) <m), 

lm<|rTJ te [ 0jL ] J 

r,2 = {z r > m (t) G K for all m < |rT] and t G [0, L]}, 
^ 3 = { max sup |(l,^ m (i))-(l,^ m (0))|<-j, 

^ 5 ^n{ m ^ TJ w,(^(-),^)<i}, 

^ 6 c(sup sup |£^(A)-a r (t- S )^(A)|<ij, 

[s,t]C[0,£] 

n; iL (M, k, c L ) = n 5 "- 1 n n r ' 2 n ^ 3 n n ^ 5 n n£ 6 , 

where for each fixed n G N, the events {f2^' 6 : r G 7?-} are chosen so that 
lim r _ >00 P r (J}£ 6 ) = 1 (see Lemmas 4.7 and 4.9). 
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Fix M, K and collections {Cl-L > 1}. Then for each L > 1, the events 
{Q^ l (M,K,Cl) -r G TZ, n G N} form an array indexed by TZ x N. The next 
lemma asserts that M, K and {Cl ■ L > 1} can be chosen so that for each 
L > 1 , the tail of every TZ- indexed column in this array has arbitrarily high 
probability. 

Assume (A) and let T > 1 and < r] < 1. 

Lemma 4.16. There exist M* > 0, a compact K* C M and collections 
{C£ : L > 1} such that for each L > 1 and each fixed n6N, 

(4.74) hminf P r ^ L (M* ,K* ,C* L )) > 1 - \. 

r — >oo ' 2 

Proof. Use Corollary 4.3 and Lemma 4.10 to choose M* and K* so 
that for all L > 1 and each fixed n, 

(4.75) limmf P r (n r ' 1 n !T' 2 n ^' 3 ) > 1 - ~. 

Let L > 1. By Lemmas 4.12, 4.14, 4.7 and 4.9, there exists C£ = ({«j},{^}) 
such that for each fixed n, 

(4.76) liminf P r (^ 4 n n ^' 6 ) > 1 - \. u 

The preceding lemma guarantees that for each L > 1, there exists a "di- 
agonal" sequence of events, indexed by r G such that the tail of this 
sequence has arbitrarily high probability. 

Assume (A) and let T > 1 and r/ £ (0,1). Using Lemma 4.16 choose a 
constant M*, a compact set K* and collections {C£:L > 1} so that (4.74) 
holds for each L and n. 



Definition 4.17. For each L > 1 and r 6 TZ, let 

HI = {n G N : P r (^, L (M*, K*, C£)) > 1 - 7?} 

and let 

(4.77) n(L, r) ={7^ A W' W 

For each L > 1 and r &TZ, define 

(4-78) n r L = n r n{L>r)jL (M*,K*,c* L ). 

Note that in the above definition, supA/"£ could equal infinity for some 
r. In this case, it is sufficient to use VL r L = JIT. L (M*,K*,C2), which is the 
reason for including the minimum with \r\ in the definition. This does not 
affect the property that n(L,r) G A/"£ for M£ /0: if m,nGM with m <n, 
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then n r n>L (M*,K*,Cl)cn r m>L {M*,K*,Cl). So, for m < n, n G N r L implies 
that m G A/"£. In particular, A/£ 7^ implies that n(L,r) G A/£. 

Also, note that f2£ may be empty for some r, possibly when A/£ = 0. 
This is of no concern; the following lemma implies that this is the case for 
at most finitely many r. 

Assume (A) and let T > 1 and rj G (0, 1). 

Lemma 4.18. For each L > 1, 

(4.79) liminf P r (ST L ) > 1 — 77. 

Proof. Fix L > 1. If n(r, L) > 1 for some r G 72, then P r (Q r L ) > 1 — 77 
by definition. Thus, it suffices to show that 

(4.80) n(r,L)—> 00 as r -> 00. 
This follows by Definition 4.17 and Lemma 4.16. □ 

Using the events fi£, define, for each L > 1, a sequence of subsets of 
D([0, oo),M) that are sample paths on [0, L) of the shifted, fluid scaled 
state descriptors {Z r ' m (-) :m = 1, . . . , L r ^"J}- Assume (A) and let T > 1 and 
77 6(0,1). 

Definition 4.19. For each L > 1 and r eTZ, let £F£ be the set of all 
C(-) G D([0, 00), M) such that for some u G f2£ and some m < [rT\ , 



(4.81) 




t G [0, L), 
t G |X,oo). 



By Lemma 4.18, Ql r L is empty for at most finitely many r. 

Assume (A) and let T > 1 and 77 G (0, 1). For each L > 1, let : r G 7£} 
be the sequence of sets defined above. Fix L > 1 and suppose that 1Z C 1Z 
is a subsequence and {C r (0 : ^ £ ^} C D([0, 00), M) is a sequence of paths 
such that C{-) 6 S£ for each f G 71.. 

Theorem 4.20. TTie sequence {C, r '(•): f G 72.} is precompact in D([0, 00), 
M). 

Proof. Let K* and C£ = ({re?}, {^fc}) be given by Lemma 4.16. By 
definition of ^£ and 0£ (O f ' 2 in particular), ( f (t) G K* U {0} for all t > 
and f G 72. Let /c G N. By (4.80), there exists such that f > implies 
that n(r,L) > fc. For R > 1, let w^(-,-) be the modulus of continuity on 
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[0,R] used in Theorem 3.6.3 of [7]. Since lim^o w j?(£(')j <5) =0 f° r each 
fixed £(•) G D([0,oo),M), there exists S' k G (0,<5 fc ) such that 

(4.82) maxw^(C f (-),^)<^ 

If f > r^, then 0/'^ 5 ~ L ^ C Sl)j 5 (see Definition 4.15). By definition of Ss r L and 
W L {$X,l(r L) ^ n P ar ticular) , we have 

(4.83) sup w L „(C ? (-),^) < sup w L „(C f (-),4) < \- 

r>r k r>r k & 

Since w^(C r (-), S' k ) is bounded above by wl_((7(-), <%) for each r <E1Z and 
A: G N and since C r (") = on [L, oo) for each f, (4.83) implies that for R> 1, 

(4.84) lim sup w' fl (C f (•), S'k) < 1™ ?up w L -(C (•). 4) = «• 

Deduce from (4.82) and (4.84) that 

(4.85) lim supw^(n-)X) = 0. 

So {C r (") : ^ £ 7£} is precompact by Theorem 3.6.3 of [7]. □ 



5. Local fluid limits. This section provides the remaining ingredients 
needed to prove Theorem 2.2. It is divided into three parts. Section 5.1 
below identifies the limit set S>l of the sequence of sets {& r L :r G 1Z} intro- 
duced in the previous section and derives several properties of the elements 
of this set. Section 5.2 establishes that the set 2$l approximates the sample 
paths of Z r,m (-) on [0, L) with high probability. Section 5.3 contains a result 
on the steady state behavior of the elements of au important ingredi- 
ent for proving state space collapse of the diffusion scaled state descriptors 
{Z r (-)}. The results of this section are combined in Section 6 to show state 
space collapse, which quickly leads to the proof of Theorem 2.2. 

5.1. Properties. Assume (A) and let T > 1 and rj G (0, 1). 

Definition 5.1. For each L > 1, let be the set of all ((■) G D([0, oo),M) 
such that there exists a subsequence IZdTZ and a sequence {C(-) : f G 1Z} C 
D([0,oo),M), satisfying C(-) G & f L for each r G H and 

(5.1) CT)^C(-) asf^oo. 

The set Qi, is nonempty by Lemma 4.18 and Theorem 4.20. Call the 
elements of local fluid limits on [0,L) of the state descriptor Z r (-); the 
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term indicates that each element £(•) £ ^x, approximates some section of a 
sample path of Z T {-) : [0, rT] — > M over a finite time interval of length L. 
The precise statement of this approximation appears in Section 5.2 below; 
the purpose of this section is to establish several properties of the elements 

Let denote the Borel subsets of H+. Assume (A) and let T > 1 and 
r?G(0,l). 

Theorem 5.2. There exist a positive constant M and a compact set 
K C M such that for each L > 1 and each £(•) £ 

(i) ((t) £ K /or 4 > 0; 

(ii) C(-) = on [L,oo); 

(iii) £(•) *s continuous on [0, L); 

(iv) (l,((t))<M for allte [0,L); 

(v) (1, (■(•)) is constant on [0,L); 

(vi) i/ z = (1, C(0)) ^ s positive, then for all t £ [0, L) and B £ ffi ' , 

(5.2) C(i)(^) = C(0)(B + (^- 1 ,t)) + a f $(B+(sz-\s))ds. 

Jo 

PROOF. Let M*, K* and C£ = ({kj}, {S n }) be the objects used in Def- 
initions 4.17 and 4.19 to define W L . Let M = M* and K = K*. Fix L > 1 
and ((■) £ 3>l- By Definitions 4.19 and 5.1, there exist a subsequence IZclZ 
and a sequence {C(-) £ 7^} such that C r (") S ^£ for each f and such that 
(5.1) holds. To ease notation, assume without loss of generality that 1Z = 1Z. 
Properties (i) and (ii) follow immediately from Definitions 4.17, 4.19 and 
5.1. Property (iii) follows from the proof of Theorem 4.20, in particular, 
from (4.83). 

Since ((■) is continuous on [0,L), (5.1) implies that for all t £ [0,L), 

(5.3) ( T \t)^U({t) asr^oo. 
Consequently, for all t £ [0,L), 

(5.4) (l,C(t))= lim (l,C(t))<A/, 

r — >oo 

by the definition of *3> r L (O r>1 in particular). This proves (iv). Property (v) 
follows from (5.3) and the definition of W L (^' 3 in particular). 

To prove (vi), suppose that z = (1,£(0)) is positive and fix t £ [0,L) and 
A £ stf . An extension to Borel sets B £ Jrf? will be established at the end 
of the proof. For each r £ 1Z, there exists u £ W L and m < [rT\ such that, 
evaluated at to, 

C(t)(A) = C(0)(A + (S r > m (t),t)) 

(5.5) 
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rE r ' m (t) 

+- r e ^rm m m 

i=rE r <™ (0)+l 

see Definition 4.19 and the dynamic equation (4.21). By (4.83) and (5.3), 
(5-6) lim ||<l,C r '(-)> - <1,C(-)>IU = 0. 

r — >oo 

By property (v), (l,((u)) = z > for all u G [0,t\. Therefore, (5.6) and the 
bounded convergence theorem imply that for all s G [0, t) , 



lim S r ' m (s) = lim / <p((l,C(u)))du 

r— *oo r—*oo Jq 

(5.7) = HixHr'du 



= sz-\ 

For s G [0,t], define S s j = (t — s)z . Since s t— > sz~ 1 is continuous on [0,t], 
the convergence in (5.7) is uniform on [0, t]. For 5 > 0, there exists rg G 1Z 
such that 

(5.8) sup \S r s '™ — S s ,t\ < <5 for all r>r$. 
se[o,t] 

The following lemma will be needed. 

Lemma 5.3. For oil i 6 j;/ and s G [0,L), 

(5.9) CW(^a)=0. 

Proof. Fix A G and observe that 8a C <9^ n for all n G N. For all 
t G [0, L) and n G N, (5.3), the Portmanteau theorem and the definition of 
@ r L (O^' 4 in particular) yield 

C(t)(9A) < at)(9T) < limmf C (t)(&?) < ~- 
Letting n — > oo proves (5.9). □ 

Continuing with the proof of Theorem 5.2(vi), we use (5.7) and the defi- 
nition of £F£ (^' 4 i n particular) to obtain 

(5.10) \ CmA +(S r > m (t),t))-e(0)(A + (tz-\t))\^0 

as r — > oo. Deduce from (5.3), (5.10), Lemma 5.3 and the Portmanteau the- 
orem that as r — > oo , 

(5.11) C(t)(A) - C(0)(A + (S r > m (t), t)) -+ C(t)(A) - C(0)(A + (tz-\t)). 
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Let D#{&f) = {D G £/:&(d D ) / 0}. Note that D#(&/) is countable be- 
cause •&(• xl) and #(R+ x •) are probability measures. Since the function 
s | — ► S^t is strictly decreasing in s, 

£>* (S) = {s £ [0, t] : A + (S M ± 2<5, t - s) € A?(0} 

is also countable. For each integer N > 1, let = t$ < t± < • • • < tjy = £ be a 
partition of [0, i] such that tj^ ^ D$(S) for all j = 1, . . . , N — 1 and such that 
maxj<7v-i(ij^ 1 — i^) — > as N — > oo. Let I r denote the second right-hand 



term in (5.5). Then 

N-1-. rE^{tf +1 ) 

r=Y,\ E i + A W m m m {t)). 

Suppose that < U\r~ x — m < t^ +l for some r > rg, some j < N — 1 and 
some i E {rE r > m (0) + 1, . . . ,rE r > rn (t)}. Then, by (5.8), 

( 5 - 12 ) %. . ,t - 5 < ^ar>-i-m t ^ %,t + 

By the definitions (2.14), (2.15), (2.5) and (2.3), 

(vnt),rr(t)) = « - c-™,^ - (*- (^ - *o)). 

So, for r > r$, (5.12) and the inequalities 1a(- — 5, •) < 1^(-, •) < 1a(- + 5, •) 
yield 

i5(<' m (ty? m (*)) > - (s tfjt + 2^,/rr- 1 - (* - tf )) ; 

l+«^(t),ZT^(t)) < U« - (5 tf+i)t - M),^" 1 - (t - tf +1 )). 
This yields, for r > r$, 

N—l -j ^ ro (f+i) 
^ > E ^ E lA« - (fy, t + 25), /[r- 1 - (t- )); 

3=0 i=r_B r . m (t JV )+l 

(5.13) 

^ ^ E - E M< - (S t ^ t - 25)Xr~ l - (t - tf +1 )). 

3=0 i=rE r < m {tf)+\ 

Using (2.16) and (2.18), we can rewrite (5.13) as 
JV-l 

r > £ (A + (S tN t + 2S,t- if)); 

3=0 J ' J+1 

(5.14) 

iv-i 



r < £ £ r t'™ t » ( A + ( s t-t -M,t- tf +1 )). 

3=0 
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By Definitions 4.17 and 4.15, (5.14) implies that, for r > r$, 

N-l 

r > E (« r (*f+i - f W(a + (% * + 25,* - *f )) - r)- 1 ), 

i=o 

AT-l 

/ r < E (« r (*f+l " *f )^ + (V +1 , t - 25,* - *f +1 )) + n(L, rn 1 ). 
i=o 

By (4.80) and the fact that tf <£ D#(S) for all j = 1, . . . , N - 1, 

JV-l 

liminf r >«E (*f+i " *f W + (%,* + 25, * - *f )), 

j=o J 

(5.15) 

AT-l 

hmsup r < a E (*f +1 " *f )^ + 0V +1 ,t " 25, * - *f+i))- 

r— >oo . „ J+ 1 

For s G [0,t] such that s ^ D$(S), the bounded convergence theorem implies 
that as N ^ oo, 

iV-i 

E l [t M >tN )(sMA + (5 t * t + 25, t - tf )) - + (5,, t + 25, t - s)), 

. — * 1 J ' 3+1' J ' J 

(5.16) 



N-l 

E l [tfitf+i) {s)V{A + (S t N +i>t -25,t- tf +1 )) - tf(,4 + - 25, t - s)). 

j=0 



Thus, the convergence in (5.16) holds for almost every s € [0, t). Let N — > oo 
in (5.15) and conclude from (5.16) and the bounded convergence theorem 
that 

liminf I r >a t #(A+ (S st + 2S,t- s))ds, 

r— too In 

(5.17) 

limsup/ r <a / #(A + (S s ,t-2S,t-s))ds. 

i — too JO 

Let 8 — > in (5.17). Since D&(£/) is countable, both integrands in (5.17) 
converge almost everywhere on [0, t] to 



&(A + (S s , t , t-s)) = 0(A + ((t - s)z-\t - s)). 
Conclude, using a change of variables, that 

lim I r = a $(A+(sz-\s))ds. 
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Combining this with (5.11) proves (5.2) for fixed t G [0,L) and A G srf . 
To extend to all B G , let be the set of B G for which (5.2) holds. 
Observe that =2/' is a A-system: H + G stf 1 because H + G .2^; if {-B n } C si 1 
satisfies B n f 1?, then G if B\ C -B2 are elements of <s/' , then B2\B\ G 

Also, observe that srf is a 7r-system: if A\,A% G then ^1 n ^2 G 
Since srf C and the cr-algebra generated by stf is equal to jri? , it follows 
that &4' = ^ by the Dynkin 7rA-theorem (see, e.g., [1]). □ 

5.2. Uniform approximation. The next lemma states that with asymp- 
totically high probability, the overlapping sections {Z r ' m (-) :m < L r ^J} of 
the process Z r {-) are uniformly approximated on [0, L) by elements of the 
set 3> L . 

Assume (A) and let T > 1 and rj G (0, 1). Let {S>l ■ L > 1} be the subsets 
of D([0, 00), M) specified by Definition 5.1. 

Lemma 5.4. For each L > 1 and each e > 0, the set 3>l satisfies 



(5.18) liminfPM max inf sup d\Z T ' m (t), C(t)] < £ > 1 - n. 

™ \™<\rT\a-)£®Lte{o,L) wswj - y- 

Proof. Fix L > 1. By Definition 4.19, Z^ m (')l[ 0)L ) (•) G ^£ for all m < 
[rTj if and only if uj G So, by Lemma 4.18, 

(5.19) liminf P r (^ r ' m (-)1[ ,L)(-) G 0£ for all m < L^TJ ) > 1 - 77. 
Thus, it suffices to show that 

(5.20) lim sup inf sup d[C{t), C(*)] = 0. 

r-+oo fr(.) 6 #r f G [0,L) 

Suppose that (5.20) does not hold. Then there exists e > 0, a subsequence 
IZcTZ and a sequence {C r (') : ? £ 72.} with C r (") G ^£ f° r each f such that 

(5.21) inf inf sup d[( f (t),((t)} > e. 
f&nC{-)^Lt£[o,L) 

By Theorem 4.20 and Definition 5.1, there exists a further subsequence 
{fj} C 1Z and a £(•) G @l such that 

(5.22) C ^(.)^ C( .) as./ --X. 

The path £(■) is continuous on [0, L) by Theorem 5.2(iii). Therefore, 

(5.23) lim sup d[C fj (*),C(t)] = 0, 

^°°te[o,L) 

contradicting (5.21). □ 
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5.3. Convergence to steady state. This section establishes that when L 
and t are sufficiently large, is arbitrarily close to a steady state measure 
in M#, uniformly for all £(•) G The proof of this assertion exploits (5.2) 
and the constant total mass function (1, £(•)). 

Assume (A) and let T > 1 and e, rj G (0, 1). For each L > 1, let @l be the 
set of local fluid limits specified by Definition 5.1. 

Theorem 5.5. There exists L* > 1 suc/i tfiai for all £(•) G ^l*, 

(5.24) sup d[C(<),A*(l,C(*))]<e. 
te[L*-i,i*) 

Proof. Let M and K be given by Theorem 5.2. Choose K > such 
that 

(5.25) sup^Q-fTM -1 , oo) x R) < e 
and such that 

(5.26) a/ $([sM- 1 ,oo)xR)ds<e; 

Jk 

such a if exists because K is compact and (x,"&) < oo by (2.21). Choose L* 
such that 

(5.27) L* - 1 > K. 

Fix C(-) G ® L * and t G [L* — 1,L*). By Theorem 5.2(v), (1,C(0) is constant 
and so 2 = (1,£(0)) satisfies 

(5.28) A^l,C(t)) = AXl,C(0))=^. 
By Theorem 5.2(iv), 

(5.29) z < M. 

Let B C 1H + be a closed Borel set. By (5.28) and the definition of the 
Prohorov metric d[-, •], it suffices to verify the two inequalities 

(5.30) at)(B)<r e (B E )+e, 

(5.31) r e (B)<((t)m + e. 
By Theorem 5.2(vi), 

(5.32) ({t)(B) = ((0){B + (tz-\t)) + a f $(B + (sz" 1 ,s)) ds. 

Jo 

Since t>L* — l>K, (5.29) and (5.25) imply that 

C(0)(J3 + (te -1 ,*)) < C(0)([KM- 1 ,oo) x R) 
< e. 
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Substitute this into (5.32); for the second term, enlarge B to B £ and enlarge 
the domain of integration to conclude that 

pOO 

C(t)(B) <e + a ®{B £ + (sz-\s)) ds = e + $ z e {B £ ). 
Jo 

This proves (5.30). To prove (5.31), use Definition 2.1 to write 

pt poo 

$ z e (B) = a ^( J B + (sz~ 1 ,s))ds + a / ${B + (sz~\ s)) ds. 
Jo Jt 

Enlarge B to B e in the first right-hand term; use t > L* — 1 > K and z < M 
in the second right-hand term to obtain 

pt poo 

$ z e {B)<a $(B £ + (sz- 1 ,s))ds + a 0([sAf _1 , oo) x R) ds. 
Jo Jk 

By (5.26), 

•01(B) < a f $(B £ + (sz-\s))ds + e. 
Jo 

Add an extra term and use Theorem 5.2(vi) to conclude that 

K(B) <((0)(B £ + (tz-\t)) + a C ®(B £ + (sz~\s))ds + e 

Jo 

= C(t)(B £ ) + e, 
which proves (5.31). □ 

6. Diffusion limit. The work of the previous sections is now combined to 
prove state space collapse, which immediately leads to the proof of Theorem 
2.2. 

Assume (A) and let T > 1. 

Theorem 6.1. Asr^oo, 

(6.1) sup d[Z r (t),A^(l,Z r (t)}}^0. 
te[o,T] 

Proof. Let e, rj € (0, 1). It suffices to show that 

(6.2) liminfP'Y sup d[Z r (t), A#(l,Z r (t))} < e) > 1 - n. 

r^oo Vte[o,T] / 

By Theorem 5.5, there exists L* > 1 such that the set 3>l* defined in Defi- 
nition 5.1 satisfies 

(6.3) sup sup d[C(t),Atf<l,C(t))]<|. 
C(-)es> L * te[L*-i,L*) 3 
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By Theorem 5.2(iv), there exists a positive M such that for all £(•) G 

(6.4) sup (1,C(*)> <M. 

te[o,L*) 



Choose a positive 5 < e/36 such that for each z G [0, M+ 1], 

(6.5) SUp (lb^+M] x R,^e) < SU P (l[*,x+9«]xR>^ +1 > < f 



(see Definition 2.1). The map A^:M + — > is uniformly continuous on 
[0,M + 1] (Lemma 4.5). Thus, there exists a positive S\ <5 such that 

(6.6) sup d[A^z, A$y] < 5. 

y,z<M+l. 
\z-y\<8i 

For each r &1Z, define events 

nj = { max inf sup d[^ m (t), C(t)] < sX 

n5 = {d[^°(o),A tf (i,^°(o))]<5 1 }, 
n r = n{ n ir 2 . 

By Lemma 5.4, (2.30) and (2.32), 

(6.7) hminf P r (^o) >!-»/• 

Fix lu € and assume henceforth in the proof that all random objects are 
evaluated at this to. By (6.7), it suffices to show that 

sup d[Z r (t),Av(l,Z r (t))]<e. 

te[o,T] 

Note that if t G [L*/r,T], then Z r (t) = Z r (rt) = Z r > m (s) for some m < L^TJ 
and some s G [L* - 1,L*). If i G [0,L*/r), then i r (t) = Z r {rt) = Z r >°(s), for 
some s G [0, L*). Thus, it suffices to show the two inequalities 

(6.8) max sup d[Z r > m (t), A* (1, Z r ' m (t))} < e, 

m<[rT\ te[L*~l,L*) 

(6.9) sup d[Z r >°(t),A#{l,Z r '°(t))} <£. 

te[o,L*) 

Fix m < \rT\ and t G [0, L*). By the definition of Q\, there exists £(•) G i^L* 
such that 

(6.10) d[Z r ' TO (t),C(t)]<«i. 

By definition of the Prohorov metric d[-, •], this implies that 

|<i,C(f)>-(i,J r ' m (t)>|<*i<i 
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and so by (6.4) and (6.6), 

(6.11) d[A#(l,((t)),A#(l,Z r ' m (t))]<5- 
By (6.10) and (6.11), we have 

d[Z^(i), A,<1, Z r > m (t))] < d[Z r ' m (t), cm 

+ d[C(t),A^(l,C(t))] 

(6.12) 

+ d[A,,(l,C(t)),A^(l,^ m (t))] 
<2«5 + d[C(t),Atf(l,C(t))]. 
If t G [L* — 1,L*), substituting (6.3) into (6.12) yields 

d[Z*' m (t), A*(l, Zr> m {t))\ <25+ e -<e, 

which proves (6.8). To prove (6.9), it suffices to show that for m = and 
i£ [0,L*), the £(•) appearing in (6.10)-(6.12) also satisfies 

(6.13) d[C(t),A*<l,C(t)>]<|. 

By Theorem 5.2(v), i s constant. Let z = (1,£(0)) so that 

(6.14) A^(l,C(t)) = A^(l,C(0))=C- 

Fix a closed Borel set B C M + . It suffices to show the two inequalities 

(6.15) C(i)(£)<^(£ £/3 ) + " 



3' 

(6.16) #*(!?) <C(*)(£ £/3 ) + f. 
By (6.14), (6.10), the definition of O r 2 and (6.11), 

d[C(0),^]<d[C(0),jT>°(0)] 

+ d[^°(0),A 1? (l,^°(0))] 

(6.17) +d[A d (l,^°(0)),^] 

<5 1 + 5 1 + 5 

<35. 

By Theorem 5.2(vi), ((■) satisfies (5.2). Apply (6.17) to the first term on 
the right of (5.2) to obtain 

(6.18) C{t){B)<^ z e {{B + {tz- 1 ,t)f 5 ) + 2,5 + a ${B + (sz" 1 ,a)) ds. 

JO 
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Note that 

r e ((B + (tz~\t)f) < tf z e (B 3S + (tz-\t)) 

(6.19) 

+ sup (l[ x ,x+36]xM^ 



I) 

x& 



Apply (6.5) and use the relation B C B 3S C B £ l 3 to conclude from (6.18) 
and (6.19) that 

(6.20) ({t)(B)<tf z JB £ / 3 + (tz-\t)) + a f ${B £ I 3 + (sz~\ a)) ds + 35+ -. 

Jo 6 

Inequality (6.15) follows from (6.20) and Definition 2.1. To show (6.16), 
write 

d z e (B) = {) z e {B+(tz- 1 ,t)) + a / d(B + {sz- 1 ,s))ds. 

Jo 

Use (6.17) and enlarge B to B 3S to obtain 

(6.21) ^(S)<C(0)(( J B + (^- 1 ,t)) 35 )+35 + a / ^(S 35 + (sz" 1 , s)) ds. 

Jo 

Note that 

((0)((B+(tz~ 1 ,t)f)<aO)(B 35 + (tz-\t))+ sup (l[^xK,C(0)) 

and so by (6.17), 

C(0)((B + (te- 1 ,*)) 35 ) < C(0)(B 35 + (tz" 1 ,*)) 

(6.22) 

+ SUp (l[ x ,x+9S]xU,K) + 3(5 - 

Apply (6.5) to deduce from (6.21) and (6.22) that 

#t(B) < ((0)(B 35 + (te -1 ,*)) + a [^{B 35 + (sz" 1 , a)) ds + 65 + -. 

Jo 6 

Conclude from Theorem 5.2(vi) that 

(6.23) r e (B)<C(t)(B 3S ) + 6S+^ 

6 

which proves (6.16). This completes the proof of (6.13), which, together 
with (6.12), proves (6.9). □ 

Proof of Theorem 2.2. By Proposition 4.2, 

(l,i r (-)) =►£*(•) asr^oo. 
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Since the lifting map A$ is continuous (Lemma 4.5), the continuous mapping 
theorem implies that 

A#(l,Z r (-)) =>£*(•) asr^oo. 

Conclude from Theorem 6.1 and the "converging together lemma" ([2], The- 
orem 3.1) that 

Z r (-)^Z*(-) asr^oo. □ 
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